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PREFACE 


In  the  fall  of  1976  the  author  was  invited  to  render  a lecture  on  Cas 
Bearings  in  a course  entitled  "Principles  of  Trlbologv"  at  ttie  University 
of  Pittsburgh.  A number  of  texts  on  the  subject  are  available  ' <6,  37, 

38,  39,  40,  and  41).  Together,  these  publications  cover  the  field  quite 
thoroughly.  However,  being  the  products  during  an  era  ol  rapid  technolog- 
ical developments,  their  perspectives  are  somewhat  diffused.  Therefore, 
this  writing  was  undertaken  in  order  to  project  the  theoretical  foundation 
of  gas  lubrication  as  a mature  technical  discipline.  It  is  offered  as  an 
overview  reference,  supplementing  the  texts  previously  cited. 

The  encouragement  by  the  Office  of  Naval  Research  to  make  this  writing 
reach  a wider  readership  is  most  welcome.  This  distribution  as  a con- 
tract report  allows  the  writer  to  express  his  gratitude  for  the  continuing 
research  support  given  to  him  by  ONR  since  1961. 

--C.H.T.P. , May  1977 
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1.0  INTRODUCTION 

Men  have  derived  beneficial  mechanical  effects  from  air  and  other  gases  in  many 
ways.  In  most  cases,  a moving  body  is  involved,  and  the  relative  kinetic  energy 
is  converted  to  create  levitation  or  to  perform  useful  work.  Viscous  shear 
would  contribute  to  parasitic  losses  directly  or  would  cause  the  generation  of 
turbulent  flows  which  in  turn  would  precipitate  dissipation.  Compressibility 
of  the  gaseous  medium  is  often  an  essential  feature  of  the  useful  mechanical 
process;  e.g.  expansion  work  in  heat  engines  and  propulsion  nozzles.  Fluid 
film  lubrication  is  an  exceptional  mechanical  process  in  that  viscous  shear 
stress  contributes  directly  to  the  useful  function  of  developing  a load  capac- 
ity. Although  viscosity  also  causes  bearing  friction,  the  equivalent  lift-to- 
drag  ratio  of  a typical  ^hydrodynamic  wedge**  is  in  the  order  of  1000  to  1, 
which  compares  favorably  to  a high  performance  wing.  Lubricant  compressibility 
is  the  distinctive  feature  of  gas  bearings  in  contrast  to  the  more  common  liquid 
lubricated  bearings.  Although  basic  concepts  such  as  the  ^hydrodynamic  wedge 
ore  still  applicable  to  gas  bearings  despite  of  lubricant  compressibility,  many 
additional  features  in  gas  bearings  are  unique  and  require  separate  attention. 

/V- 

1 

The  potentials  for  large  scale  industrial  application  of  gas  bearings  were 
recognized  in  the  late  1950' s.  Advocates  of  gas  lubrication  have  emphasized 

• The  gaseous  lubricant  is  chemically  stable  over  a wide  temperature 
range. 

• Atmospheric  contamination  is  avoided  by  gas  bearings. 

• The  viscosity  of  a gas  increases  with  temperature  so  that  the 
heating  effect  in  overloading  a gas  bearing  tends  to  increase 
the  restoring  force  to  overcome  the  overload. 

• A gas  bearing  is  more  suitable  for  high  speed  operation. 
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• There  is  no  fire  hazard. 

• Use  of  gas  bearings  can  reduce  thermal  gradient  in  the  rotor 
and  enhance  its  mechanical  integrity, 

• For  high-speed  applications,  the  gas  bearing  is  Inherently  more 
noise  free  than  the  rolling-element  bearing. 

• System  simplicity  is  enhanced  by  use  of  self-acting  gas  bearings 
which  do  not  require  cooling/circulation  equipment. 

These  optimistic  views  must  be  tempered  with  more  subtle  engineering  con- 
siderations before  one  can  confidently  substitute  gas  bearings  for  more 
conventional  oil-lubricated  bearings. 

Surface  Chemistry  - Modern  lubricating  oils  contain  additives  of  long  chain 
molecular  species  which  form  low  friction  surfac:  coatings  on  bearing  sur- 
faces. Consequently,  sliding  friction  in  oil  bearings  is  minimal  even  when 
full  separation  of  the  bearing  surfaces  by  the  fluid  film  is  not  realized. 
Boundary  lubrication  of  gas  bearings  is  essentially  not  possible*. 

_3 

Lift-off  Speed  - Viscosity  of  gases  is  less  than  10  times  that  of  lubri- 
cating oils.  Thus,  for  a self-acting  gas  bearing,  the  lift-off  speed  is 
proportionally  higher  than  its  oil-lubricated  counterpart.  Friction  heating 
at  the  sliding  contact  prior  to  lift-off  is  accordingly  more  intense  and  is 
therefore  more  likely  to  cause  material  damage  in  self-acting  gas  bearings. 

Chemical  Deterioration  - The  gaseous  medium  is  chemically  quite  stable  at 
ordinary  environments.  Lubricating  oils,  in  contrast,  would  decompose  at 


Boundary  lubrication  is  sometimes  used  in  gas  bearing  gyroscope  to  ensure 
minimal  starting  friction  for  prolonged  usage.  The  friction  coefficient 
of  such  surfaces  is  typically  > 0.30f^J  while  that  of  conventional  bearing 
surfaces  in  the  presence  of  oil  is  usually  < O.lOt^). 


F 
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a relatively  low  temperature.  For  most  oils,  coking  would  take  place  below 
300“C.  In  the  event  of  Inadvertent  malfunction,  sliding  friction  can  cause 
considerable  temperature  rise.  If  the  coking  temperature  is  reached,  the 
endothermic  chemical  process  would  function  as  a sink  to  moderate  further 
temperature  rise  and  thus  would  prevent  material  damage  of  the  bearing 
surfaces. 

These  three  factors  combine  to  make  gas  bearings  more  susceptible  to 
mechanical  damage  and  thus  preclude  widespread  application  of  gas  bearings 
in  heavy  duty  equipment.  The  same  considerations  also  exert  a dominating 
influence  in  the  choice  of  satisfactory  materials  for  gas  bearings. 
Beneficial  use  of  gas  bearings  must  be  predicated  on  the  avoidance  of 
these  limiting  factors. 

Nomenclature  for  mathematical  symbols  are  separately  summarized  at  the  end 
of  each  chapter.  Where  different  definitions  for  the  same  symbol  are  used 
in  different  sections,  the  section  numbers  are  indicated. 


r 
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2.0  MATHEMATICAL  THEORY  OF  GAS  LUBRICATION 

Theory  of  gas  lubrication  is  generally  regarded  as  an  extension  of  the  fluid 

[3  ] 

film  lubrication  theory  of  Reynolds,  which  was  originally  formulated  for  an 
incompressible  lubricant.  TTie  crucial  issue  is  concerned  wi th  an  appropriate 
account  of  the  density  variation  within  the  lubricant  film  such  that  the 
basic  principles  of  thermodynamics  are  satisfied  to  a degree  consistent  with 
the  approximations  already  invoked  in  momentum  considerations. 

In  certain  ways,  analysis  of  a gas  bearing  is  simpler  than  that  required  for 
liquid  bearings.  As  to  be  demonstrated  below,  the  temperature  in  a gas  bear- 

14] 

ing  film  may  be  regarded  to  be  a constant  , even  though  viscous  heating 
necessarily  cause  some  temperature  rise  above  that  of  the  bearing  surfaces. 
Since  the  viscosity  coefficient  of  most  gases  is  dependent  solely  on  tempera- 
ture, an  isovlscous  approximation  is  quite  satisfactory  for  studying  gas 
bearings.  In  a liquid  bearing  film,  the  isovlscous  approximation  is  less 
reliable.  The  gas  bearing  film  is  inherently  a single  phase  constituent. 
Irrespective  of  the  local  pressure  level  relative  to  the  ambient,  the  gaseous 
lubricating  film  remains  a homogeneous  medium.  However,  in  a liquid  film 
bearing,  it  has  been  an  established  empirical  knowledge  that  a homogeneous 
liquid  state  is  assured  only  when  the  local  pressure  is  near  or  above  the 
atmospheric  pressure^^^.  I4iere  the  pressure  tends  to  become  subambient  in 
a self-acting  liquid  bearing  film,  a two-phase  flow  structure  is  prevalent. 

In  fact  a totally  rigorous  treatment  of  this  aspect  of  the  liquid  lubricant 
film  is  yet  to  be  demonstrated. ^ 

On  the  other  hand,  the  kinematic  boundary  condition  at  bearing  surfaces  need 

more  attention  in  gas  bearing  analysis.  It  is  not  unusual  for  a gas  bearing 

film  thickness  to  be  in  the  order  of  1 um  or  less.  Under  such  a condition,  a 

velocity  slip  at  the  bearing  surface  may  not  be  negligible ^ , especially  if 

the  atmosphere  is  rarefied.  Furthermore,  the  common  assumption  of  imperme- 

[8  ] 

ability  of  the  bearing  surface  may  not  be  valid  . This  is  particularly 
true  for  sintered  materials. 
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Still  one  other  complication  in  gas  lubrication  theory  arises  at  cryogenic 

temperature  levels.  The  local  Mach  number  may  be  quite  large,  such  that 

phenomena  akin  to  those  in  high-speed  gas  dynamics  involving  fluid  inertia  may 
[9] 

become  evident  . In  the  following,  the  classical  isothermal  gas  lubrication 
theory  will  be  developed  with  the  continuum  point  of  view. 

2.1  Scaling  Rules  in  Thin  Film  Viscous  Flow 

The  fluid  dynamic  point  of  view  of  gas  lubrication  concerns  a film-like 
domain  which  may  be  described  by  a small  distance  measure  from  a smooth 
reference  surface.  The  fluid  in  question  is  physically  characterized 
by  the  thermodynamic  properties  of  gaseous  media. 

The  geometrical  properties  of  the  domain  are  characterized  by  two 
scales.  The  small  scale,  C,  is  the  measure  of  the  bearing  gap.  The 
large  scale,  R,  is  the  global  measure  of  the  surface,  including  its 
overall  extent  as  well  as  its  radius  of  curvature.  The  ratio 

e = C/R  <<<  1 (2.1) 

indicates  thinness  of  the  film. 

Let  V be  the  scale  of  the  velocity  field  in  concern.  If  n is  a local 
unit  normal  vector  of  the  reference  surface,  and  if  V is  a local 
velocity  vector  of  the  flow  field,  then  lubrication  fluid  dynamics 
is  further  characterized  by  the  condition. 

|V  • n|  = 0{eV}  (2.2) 

U = V - V • n (2.3) 

is  then  the  projection  of  the  velocity  vector  onto  the  tangent  plane 
of  the  reference  surface.  These  geometrical  concepts  are  Illustrated 
in  Figure  2.1. 


Fig.  2.1 

Thin  Film  Between  Nearly  Parallel  Surfaces 
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2.2  Momentum  Conservation 

If  a local  curvallnear  coordinate  system  Is  used  such  that  the  gradient  operator 
becomes 


V7  V57  4.  n ^ (2.4) 

V ' V + " -57 

where  y Is  the  distance  measured  from  the  reference  surface;  Is  the  two- 
dimensional  surface  constrained  gradient  operator;  then  the  momentum  equllitoriura 
conditions  for  either  liquid  or  gaseous  lubrlc.ition  flows  are 


V 


p + 


3 f 


iU)  = 


0 { 


PV 


£l.  yv  , 


(2.5) 


I 1^  I = 0 {-^}  (2.6) 

cR 

T is  a time  scale  representative  of  the  duration  in  which  fluctuations  of  the 
flow  field  take  place.  The  right  hand  side  of  e.q.  (2.5)  may  be  discarded  pro- 
vided one  is  willing  to  neglect  the  numerical  values  of 


,2(£V_R) 


(2.7) 


then 


3 , 3U. 

37  ^"37^ 


(2.8) 


Now  eq.  (2.6)  shows  p may  vary  by  0{-^}  along  n whereas  it  may  vary  by 

R € R 

along  the  reference  surface,  then,  since  0{e^}  is  already  omitted  from  eq.  (2.5) 


* Note  that  one  may  allow  e 


pVR 


0{e  with  0<a<l  so  long  as  ^ is  small. 


This  is  the  basis  for  the  usefulness  of  the  turbulent  lubrication  theory 

2 

[10,11].  An  auxilliary  requirement  for  eq.(2.8)  to  be  valid  is  |£^|  < < 1; 


this  is  peculiar  to  a gaseous  lubricant 


[9]. 
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In  writing  eq.  (2.8),  one  is  also  entitled  to  assert 


^ - 0 
3y 


(2.9) 


If  the  film  thickness  is  bounded  by  y = (h^,  h2)  and  (U^,  the  boundary 

conditions  of  U,  then  eq.  (2.8)  can  be  integrated  twice  to  yield 


Il(y)  ^ ^ Ij(y)l2(h2) 

u (y)  = u,  + -t-tttt  - u,}  - {-~T  /u  " i2(y)>'f  p 


■'1  " I^(h2)  ^^2 


Ij(h2) 


where 


■X 


y IL 


(m-l) 


la. 


; m = 1,  2 


(2.10) 


(2.11) 


n being  the  dummy  coordinate  along  n.  Thus  far  the  question  of  compressibility 
has  not  been  brought  up.  Therefore  eqs.  (2.10)  and  (2.11)  are  equally  applicable 
to  liquid  and  gaseous  fluids. 

2.3  Mass  Conservation 

The  continuum  mass  conservation  condition  can  be  expressed  over  the  full  film 
thickness  as 


(ph)  + 2 * (pUh)  + Zi)  = 0 


where 


— r “T  r 

ph  =1  pdy;  pUh  =J  pUdy 


(2.12) 


(2.13) 


and  is  the  total  mass  efflux  through  both  bearing  surfaces.  Eq.  (2.12)  may 
be  formally  derived  by  direct  integration  of  the  differential  continuity  equa- 
tion with  respect  to  y.  In  so  doing,  one  must  pay  attention  to  the  functional 
dependence  of  (hj^,  h2)  on  both  t and  the  surface  coordinates  and  accordingly 
impose  Leibnitz's  rule  for  differentiating  integrals  which  have  variable  limits. 
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2.4  Energy  Conservation 

With  a liquid  lubricant,  p is  constant,  and  since  it  appears  homogeneously  in 
eq.  (2.12),  it  can  be  dropped.  Whereupon  further  integration  of  eq.  (2.11) 
with  respect  to  y (assuming  a constant  viscosity)  and  subsequent  substitution 
into  eq.  (2.12),  without  p,  would  yield  the  classical  incompressible  Reynolds 
equation. 

For  a gas,  density  depends  on  both  pressure  and  temperatlon  according  to  the 
equation  of  state: 

P = f(p,  T)  (2.14) 

The  thin  film  condition  has  justified  the  approximation  of  y-independence  of  p, 

[41 

thus  the  question  of  y-dependency  of  p is  reduced  to  that  of  T . 

For  most  lubrication  problems,  the  perfect  gas  law  is  suitable,  so  that 

P = ^ (2.15) 

where  is  the  gas  constant.  Temperature  in  the  lubrication  film  must  satisfy 
the  energy  equation: 

pC  (^  + V • V)T  - V . (icVT)  = + V • V)p  + (2.16) 

V 3t  p dt  dy  dy 

The  two  terms  on  the  left  of  eq.  (2.16)  are  respectively  convective  (and  transi- 
ent) and  conduction  cooling  effects  respectively,  while  the  two  terms  on  the 
right  hand  side  are  respectively  flow  work  and  viscous  dissipation  effects. 

(Note  that  the  thin  film  approximation  has  allowed  considerable  simplification 
of  the  viscous  dissipation  term). 

Assuming  that  the  thermal  gradient  in  the  bearing-journal  structure  is  negligible, 
then  the  extent  of  temperature  variation  in  the  lubricant  film  depends  on  the 
balance  between  the  cooling  terms  on  one  side  and  the  flow  work  and  viscous  dis- 
sipation terms  on  the  other  side.  The  flow  work  term  is  seen  to  be  of  the  order. 


' 
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upon  estimating  pressure  variation  according  to  eq.  (2.10) 

uv  . 


(2.17) 


which  is  similar  to  that  of  the  dissipation  term.  The  relative  impotance  of 
transient  thermal  capacity  to  conduction  is 


pC  2d2  ^ 20’ 

(_V)  ^ (-I^)(/)  (^- 

T < K C UT 

P 


(2.18) 


pc 


(2.7),  one  can  thus  justify  neglecting  the  transient  thermal  capacity  term. 
Similarly,  the  relative  importance  of  convective  cooling  is 


PC  V 2d  2 ^ 2 i.D 

(.dRl)  = 

^ R ^ K ^ X U 

p 


(2.19) 


which  is  again  negligible'.  Thus,  the  temperature  field  in  the  lubricant  film 
is  determined  by  conduction  cooling  of  the  heat  arising  from  viscous  dissipation 
and  flow  work.  The  fractional  temperature  rise  within  the  lubricant  film  can 
thus  be  estimated  as 


PC 


V' 


2 pc 

, 

Cl  < 

P 


(y-l)M- 


(2.20) 


Consequently,  if  the  "bearing  Mach  number"  is  small,  then  temperature  variation 
across  the  lubricant  film  thickness  can  be  neglected.  Accordingly,  one  commonly 
accepts  the  isothermal  approximation  in  the  analysis  of  gas  bearings. 


2.5  Isothermal  Gas  Bearing  Theory 

Substituting  eq.  (2.15)  into  eq.  (2.12)  and  dividing  out  the  constant  JJT,  one 
obtains 

~ [p(h2-h^)]  + -2'  • [p  r ^ Udy]  ^ = 0 

‘’1 


'.M  > 


(2.21) 
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Now,  since  gas  viscosity  depends  only  on  temperature,  it  can  also  be  regarded  as 
constant.  Therefore  eq.  (2.11)  becomes 


i^(y)  = 


(y-h^) 


i2(y) 


(y2-h^?) 

~2v 


then  eq.  (2.10)  becomes 


Vp 


/ 1 \ /IT  tt 


V p 


(2.22) 


where  h = h2  - h^^.  The  two  groups  I'nclosed  by  brackets  are  respectively  the 
Couette  (linear  profile)  and  the  Poiseuille  (parabolic  profile)  components  of 
the  velocity  field. 

Integrating  across  the  film  thickness: 

y 

-y , 1 ->  >■  P 3 

Udy  = - (U^  + U2)(h)  - ^ li  (2.23) 


(h^p'^p)  + -||  = 0 (2.24) 


substitute  into  eq.  (2.21)  one  finally  finds 


^ (pi.) 


[|(U,  y 02>pi.]  - 15^  yy 


2.6  Coupling  Effects 

While  Equation  (2.24),  together  with  appropriate  boundary  and  initial  condi- 
tions, completely  describes  the  mathematical  conditions  to  be  satisfied  by 
the  film  pressure  in  a gas  bearing,  it  is  often  not  possible  to  complete 
its  solution  without  considering  coupling  effects  which  require  simultaneous 
solution  with  another  system  of  equations.  These  coupling  effects  most 
commonly  appear  in  one  or  more  of  the  following  three  forms: 
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Surface  flux,  is  a function  of  the  film  pressure. 

The  bearing  film  thickness  depends  on  the  film  pressure 
due  to  either  flexibility  (foil  bearings)  or  compliance 
of  the  bearing  surface. 

The  bearing  film  thic’.cness  depends  on  the  kinematics  of 
bearing  components  which  depends  on  the  global  effects 
of  the  film  pressure. 


The  nonlinear  appearance  of  (p,  h,  <J>)  in  Equation  (2.24)  makes  the  overall 
mathematical  problem  of  the  total  system  quite  unwieldy.  Fortunately,  the 
"external"  problem  is  usually  linear.  It  is  often  possible  to  find  the 
"Green's  function"  of  the  external  problem  so  that  such  a coupling  effect 
can  thus  be  expressed  as  a linear  integral  equation.  Thus,  for  a gas 
bearing  with  a porous  wall,  the  Green's  function  for  the  isothermal  viscous 
gas  flow  through  a porous  medium  can  be  found  so  that  one  can  write 


1 


(2.25) 


G^  is  the  corresponding  Green's  function  and  dA  is  an  area  element  of  the 
porous  surface.  Similarly,  in  the  case  of  a flexible  or  compliant  bearing 
surface,  a suitable  elasticity  equation  can  be  first  solved  to  yield  a 
Green's  function;  then  the  elastic  deformation  of  the  compliant  surface  can 
be  written  as 

Ah  = ' p G dA  (2.26) 

h 

The  Green's  function  problem  in  either  case  is  a major  topic  by  itself, 
for  the  present  purpose  it  suffices  to  indicate  that  either  Equation  (2.25) 
or  Equation  (2.26)  would  be  solved  together  with  the  Isothermal  gas 
lubrication  equation. 

When  the  kinematics  of  bearing  components  is  involved,  time  dependence 
becomes  an  inherent  feature  in  addition  to  coupling.  In  many  cases,  time 
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dependence  may  be  treated  by  a pertarbation  of  li;  i.e. 


h = h (r)  + 6h(r , t) ; |6h|  <<  |b  I 

o ’ ' ' o 


(2.27) 


r being  the  local  position  vector.  Accordingly,  the  film  pressure  can  be 
assunu'd  to  be 


= Po(r)  + 6p(r,  t);  I.SpI  <■<  Ip^J 


(2.28) 


Then  liquation  (2.2A)  for  impermeable  surfaces  can  be  separated  into  a steady- 
state  problem: 

3 

= 0 (2.29) 


h 


2 ' • (U  + I'  ) p h - P P 
2 1 2 o o 1 2u 


and  a dynamically  perturbed  problem: 


1 


li 


3 

o ->■ 


^o  ^ ■ ^2  ^ ^^2^^'’  ■ rk 


= -Pof"-2  • U|(l5,  .V'’o-T^-^oi  Pol*»- 


(2.30) 


liquation  (2.29)  must  be  solved  first;  its  nonline.ir  character  requires  an 
iterative  process  in  the  numerical  procedure.  T,qn..tion  (2.30)  is  linear 
with  respect  to  the  dependent  variable  Ap.  5h  is  treated  as  the  driving 
function,  while  the  steady-state  solution  p^^  appears  in  various  coefficients. 
The  dynamically  perturbed  problem  will  he  treated  more  thoroughly  in  Section 
6. 


2.7  Global  Bearing  Characteristics 

In  performing  a gas  bearing  analysis,  one  is  ultimately  concerned  with  the 
computation  of  flow  rate  as  well  as  the  force  and  moment  vectors  exerted  bv 
the  bearing  film.  The  flow  rate  is  mainly  a concern  of  externa  11 y- 
pressvirized  bearings  in  order  to  establish  the  flow  supplv  requirement. 

The  flow  calculation  requires  knowledge  of  the  film  inlet  pressure,  which 
is  obtained  by  establishing  continuity  netween  the  supply  restrictor  and  the 
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film  inlet  flux.  Force  and  moment  vectors  are  required  not  only  for  t lie 
nominal  operating  condition  but  also  for  perturbation  effects  to  ensure 
static  and  dynamic  stability. 

At  the  nominal  operating  condition,  static  equilibrium  is  maintained  between 

the  external  load  and  the  Integrated  film  force;  and,  in  the  case  of  a 

self-acting  bearing,  also  between  the  drive  torque  and  the  viscous  friction 

in  the  bearing.  Both  requires  consideration  of  the  hydrodynamic  stresses 

at  the  bearing  surface.  Due  to  the  thin  film  character  of  lubrication 

problems,  the  shear  stress  is  several  orders  of  magnitude  smaller  than  the 

film  pressure;  the  scale  for  the  shear  stress  being  qV/C  while  that  for  the 

2 

film  pressure  being  yVR/C  . Therefore,  in  the  computation  of  tlu>  bearing 
force  vector,  the  shear  stress  is  commonly  neglected.  However,  in  the  com- 
putation of  the  bearing  moment  vector — e.g.,  bearing  friction  torque — 
sometimes  the  film  pressure  is  inherently  of  no  consequence;  therefore,  the 
contribution  from  shear  must  be  retained.  Thus,  a distinction  must  be 
made  between  a pressure-moment  (e.g.,  pitch-moment  of  a til t ing-p.id  hearing) 
and  a shear-torque. 


The  pressure  force  acting  on  an  element  of  bearing  surface  dA  is 


(2.  n ) 


n is  the  outward  normal  unit  vector  of  the  surface  element . 
acting  on  the  surface  element  is 


P (U--U, ) 


Y Z p]dA 


The  shear  K'rce 


(2.32) 


To  allow  for  the  contribution  of  the  ambient  back  pressure,  the  bearing 
force  vector  is  thus 


^ = ; dF  - ; np  dA 

J P J a 


n(p-p  )dA 

a 


(2.33) 


r 


1 
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The  pressure-moment  Is 


M = , rxn(p-p  )dA 

P J J 3 


(2.U) 


r is  the  position  vector  of  the  surface  element  measured  from  the  origin  of 
the  chosen  coordinate  system.  The  shear  torque  is 


f r 


h 


^h  ” '2  Z'P}^^ 


(2.35) 


If  one  is  interested  in  the  static  equilibrium  condition,  then  p would  he 

o 

used  for  p in  the  above  equations. 


To  calculate  perturbation  effects  on  the  force  vector  and  the  pressure 
moment,  6p  would  be  used  instead  of  (p-p  ) in  Kquat  ions  (2.31)  and  (2.34). 

3 


L. 


/ '■ 
- 


-*d  - > 
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2.8  Nomenclature 


Symbol 

C 

C 

P 

C 

V 

-► 

F 

-> 

F 


Im 

M 


r 

R 

Si 

t 


Def inition 


Representative  scale  for  film  thickness 
Specific  heat  at  constant  pressure 

Specific  heat  at  constant  volume 

Fluid  film  force  vector 

Force  vector  due  to  fluid  film  pressure 

Force  vector  due  to  wall  shear 

Green's  function  for  surface  deformation,  Kq . (2.26) 

Green's  function  for  surface  flux  at  porous  wall, 

Eq.  (2.25) 


Bearing  gap,  h2"^’i 


Elastic  deformation  of  compliant  surface 
Normal  coordinate  integrals,  Eq . (2.11) 
Mach  number 

Moment  vector  due  to  fluid  film  pressure 

Outward  unit  normal  vector  of  bearing 
Fluid  pressure 

Surface  position  vector 


V 


r 


Del'  initlon 


Velocity  Vector 


Coordinate 


Normal  coordinate  measured  from  reference  surface 


Numerical  coefficient 


Ratio  of  specific  heats 


Perturbation  of  ( ) from  steady  state 
Three-dimensional  spatial  gradient  operator 


Subscripts 


Two-dimensional  surface  spatial  gradient  operator 

Film  thickness  parameter,  C/R 

Dummy  variable  of  y for  integration 

Thermal  conductivity 

Viscosity  coefficient 

Density  of  gas 

Representative  scale  for  time 
Bearing  wall  mass  flux 


Pertaining  to  the  steady  state 
Referring  to  lower  bearing  surface 
Referring  to  upper  bearing  surface 
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3.0  GOVERNING  DIFFERENTIAL  EQUATIONS  FOR  COMMON  BEARING  CONFIGURATIONS 
The  most  common  bearing  geometries  include  those  of 

• slider  bearing 

• journal  bearing 

• thrust  bearing 

• conical  bearing 

• spherical  bearing 


For  each  of  these  geometries,  the  operator  assumes  a special  form.  In  the 

- ► 

following,  the  terms  associated  with  ^ in  eq.  (2.24)  are  developed  for  each 
of  these  five  bearing  geometries. 

3.1  Slider  Bearing  (Figure  3.1) 

The  simplest  slider  bearing  is  one  with  a rectangular  plan  form  and  riding  on  a 
flat  moving  surface  which  has  a uniform  rectilinear  velocity  which  is  directed 
along  one  side  of  the  slider.  A very  important  application  of  the  gas  lubri- 
cated slider  is  the  read/write  head  of  a disc  or  drum  type  computer  memorv.  The 
geometry  of  a rectangular  head  riding  on  a rotating  disc  is  more  general,  this 
will  be  used  as  the  model  in  the  following  derivation. 


— — - ' 

1 

r 

[/, 

[1 f 

Fig.  1.1 


Slider  on  a Rotating  II  i sc 


As  shown  in  Figure  3.1,  one  may  describe  the  present  geometry  with  a set  of 
Cartesian  coordinates  (x,  z)  such  that  the  plan  form  of  the  slider  is  given  by 


f .<  X < f ; 


b b 

- 2 ^ z ^ 2 


f TJr'  ' -*4  - - 


-I  Si- 


Let  (1,  k)  be  the  unit  vectors  of  this  coordinate  system.  l.ie  head  itself  may 
have  a sliding  velocity  corresponding  to  an  actuation  operation: 


head 


= i u-  + kw„ 


0.2) 


Let  the  disc  center  be  located  at  (X,  Z)  in  the  chosen  coordinate  system,  then 
the  disc  velocity  is 


Uj.  = iu(Z-z)  + k(i)(x-X) 
disc 


(3.3) 


Since  the  coordinate  system  is  fixed  to  the  center  of  the  plan  form  of  tlie 
slider,  then  considering  the  upper  and  lower  surfaces  of  the  bearing  film  to  he 
respectively  the  head  and  the  disc,  one  has 


*^1  " ^disc  ' ^\iead  " k[u3 (x-X)-W2]  ; U2  = 0 


(3.4) 


When  the  head  is  at  a large  distance  from  the  disc  center,  (x,z)  may  he 
neglected  in  comparison  to  (X,  Z)  in  eq.  (3.5).  That  is,  the  disc  velocity 
may  be  approximated  by  the  tangential  disc  velocity  at  the  head  centir. 


3.2  Journal  Bearing  (Figure  3.2) 

The  journal  bearing  is  characterized  by  a cylindrical  surface.  The  curva- 
ture of  the  surface  has  significance  mainly  in  the  film  shape  of  the  eccentric 
gap  and  in  projecting  the  pressure  field  to  form  the  resultant  bearing  force. 
From  the  standpoint  of  the  differential  equation,  the  cylindrical  surface  can 
be  simply  described  by  a Cartesian  coordinate  system  (0,  z)  such  that 
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^9  R30  ^ 3z 


(3.7) 


Each  surface  may  rotate  about,  as  well  as  translate  alcjng,  the  journal  axis. 


(3.8) 


R being  the  journal  radius.  Accordingly 
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V ■ S2)pli]  = (uJj  + 1^2^30^ 


(3. 10) 


3.3  Thrust  Bearing  (Figure  3.3) 

The  overall  configuration  of  the  thrust  bearing  is  in  the  most  part  similar  to 
that  of  a slider.  However,  the  plan  form  may  be  either  an  annulus  or  a sector 
thereof,  so  that  the  natural  coordinate  system  is  that  of  the  cylindrical  polar 
coordinates,  with  the  origin  coincident  with  the  axis  of  rotation.  Writing 
(i^,  ig)  for  the  unit  vectors  along  the  radial  and  circumferential  directions. 


Ui  = igtw^;  U,  = i0ru)2 


(3.11) 


V = ir  ^ ^ 


(3. 12) 


then  r is  the  local  radial  coordinate.  Due  to  the  inherent  curvature  of  the 
coordinate  system,  for  any 


(3.13) 


i _ 1 ^ !fe 

r ; 9r  39 


(3.  lA) 


it  ' 


tl' 

**4 


-2'J- 


Consequently , 


1/  a.  , 3 (Ph) 

2(“l  2^  30 


(3.15) 


r3. 16) 


3.4  Conical  Bearing  (Figure  3.4) 

Geometrically,  the  conical  bearing  is  naturally  described  by  the  conical  polar 
coordinates  (R,  ij)).  However,  for  a description  of  the  sliding  speed,  as  well  as 
overall  mechanical  effects,  the  cylindrical  polar  coordinates  (r,9)  are  more 
useful.  These  two  coordinate  systems  are  related  to  each  other  through  the  lialf 
apex  angle  7 : 

(3.17) 


r = RslnY; 


= 9sinY 


Except  for  the  expression  of  (U^’  ^2),  equations  (3.15)  and  3.16)  would  be  an- 
plicable  upon  replacing  (r,  0)  by  (R,  <p) . The  sliding  velocities  are,  however 


U,  = 1 , Rsinyw, ; 
19  i 


U.,  = i Rsinyco,, 
^ (h  z 


(3.18) 


Thus 


'7^'  [■|(U^  + ^ 2 ^“1 


3(ph)  1 . . 3(ph)  I 

“ 2 ^“1  '"2^  “1?“  I 


(3.19) 


I V-  (h3p\^p)  = 


R If)  + ^‘’'P 


^ [sinV  ^ (rh^p  |£)  F ^ (h3p  |£)] 


(3.20) 


3.5  Spherical  Bearing  (Figure  3.5) 

The  spherical  surface  is  commonly  described  in  terms  of  ((f,  0)  which  are 
respectively  the  co-latitudinal  and  the  longitudinal  angular  coordinates.  The 


L 


■ We/, 


A 


A 


direction  of  the  polar  axis,  k,  would  be  oriented  to  preserve  symmetry  in  the 
bearing  geometry;  i.e.,  a bearing  border  is  described  by  a small  circle  ((|) 
const.).  The  triad  of  unit  vectors  of  the  chosen  coordinate  system,  (i^,  i„  , 
n) , is  defined  in  terms  of  ((J),  0)  with  respect  to  the  reference  Cartesian  triad 

(t,  j,  k) . The  spherical  bearing  is  often  used  to  allow  "free  orientation"  of 
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the  "spin  vector".  Thus  the  spin  vector,  oj,  may  not  coincide  with  the  polar 
axis,  k.  In  general,  the  orientation  of  u>  is  given  by  the  angular  parameters 
(a,  6)  with  respect  to  the  Cartesian  reference.  These  relationships  are  il- 
lustrated in  Figure  3.5.  The  spin  vector,  and  the  associated  surface  velocity 


, i [sinacosifcosCB-O)  - cosasin(}>] 


+ n [sinaslnifcos(6-6)  + cosacosif] 


if  = u)R' 


i^  slnasin(3-0) 

) ^ 

i.  + i.  [cosasinf)  - sinacos(fcos(6-6)] 

u 


The  surface  gradient  operator  is 


and  divergence  of  f = + igfg  is 


V - ^ ^ (sin<ff^)  + 


(3.22) 


i 3. 23) 


(.3.24) 


Accordingly 


(3.25 


(3.26) 
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3.6  Useful  Coordinate  Transfcrmatlons 

->•  ^ 

It  is  seen  above,  that  due  to  coordinate  curvature,  the  expressions  of  ”2  • (h^p  p) 
respectively  for  the  thrust  bearing,  the  conical  bearing,  and  the  spherical  bear- 
ing become  quite  cumbersome.  It  is  sometimes  useful  to  perform  a coordinate 
transformation  to  facilitate  analysis. 


For  thrust  and  conical  bearings,  the  required  transformation  is 


Or,  

C = ln(r/r^);  r = r^exp{;} 


(3.27) 


where  r is  a reference  radius, 
o 


For  the  spherical  bearing,  the  required  transformation  is 


sin(}i 


Or, 


y 

i;  = In  tan  -j;  slnij)  = sech^;  cosif  = -tanc 


(3.28) 


/ ‘ 
■>* 

r 


r-  ' ■-♦4  - 1 


I 
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3.7  Nomenclature 


Symbol 

a 


b 

C 

e 


f 


F 

X 


,F 


y 


F 

r 


,F 


t 


h 


H 


Chapter/ 

Def inltlon  Section 

Slider  width  in  x-dlrection 
Slider  width  in  z-direction 
Radial  clearance 
F.ccentriclty 

Radial,  circumferential,  and  col att i tudlna 1 components 
of  f 

An  arbitrary  vector 

X-  and  y-components  of  journal  bearing  force 

Radial  and  tangential  components  of  journal  bearing  force 

Bearing  gap,  h^-h^ 

h/C 


Unit  vector  along  x-axis 

Radial  and  circumferential  unit  vectors 


j 


Unit  vector  along  y-axis 


-> 

k 

P 

r 

r 

o 

R 

R 

R 

u 


Unit  vector  along  z-axls 
Fluid  pressure 

Cylindrical  radial  coordinate 
Reference  radial  coordinate 

Journal  radius  3.2 

Conical  radial  coordinate  3.4 

Spherical  radius  3.5 

x-component  of  velocity 

Projection  of  velocity  vector  on  reference  surface 


t 

i 

i 


I 


f 3 


r 
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Symbol 


head 


X 

X 


z 

Z 

ot 

a 

B 

Y 

-V 

€ 

0 


CO 


Chapter/ 

DeF Inltion  Section 

Traversing  velocity  of  read/write  head 

z-component  of  U 

Coordinate 

x-coordinate  of  disk  center 
Coordinate 

z-coordinate  of  disk  center 


Attitude  angle  3.2 

Inclination  of  spin  vector  from  polar  axis  3.5 

- ► 

Dihedral  angle  of  planes  containing  w and  k 3.5 

Half  apex  angle  of  cone  3.4 


Two-dimensional  surface  spatial  gradient  operator 


Eccentricity  ratio  of  journal,  e/C  3.2 

Transformed  merldianal  coordinate 
Polar  angular  coordinate 

Angular  coordinate  on  conical  surface  3.4 

Colattitudinal  angle  of  spherical  bearing  3.5 

Rotational  speed 

Spin  vector  of  spherical  bearing  3.5 


Subscripts 

1 

2 

(■) 


Referring  to  lower  bearing  surface 
Referring  to  upper  bearing  surface 
Time-derivative  of  ( ) 


r ' ‘ ^ ' 


I 


f 
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4.0  REPRESENTATIVE  SOLUTIONS  OF  SELF-ACTING  CAS  BEARINGS 

In  a self-acting  gas  bearing,  compressibility  is  manifested  mainly  in  terms  of 
the  relative  magnitude  between  the  viscous  wedge  action  and  tlie  ambient  pres- 
sure level.  At  a highly  compressible  condition,  the  film  pressure  tends  to  be 
so  adjusted  to  preserve  a Couette  character  in  the  velocity  profile  at  all  in- 
terior points  with  the  exception  of  sliding  exit  edges,  where  through  a sharp 
pressure  gradient  the  ambient  boundary  condition  is  satisfied.  Examples  illus- 
trating this  feature  are  given  below. 

4.1  Infinitely  Long  Sliders  (Figure  4.1) 

This  is  a special  case  of  3.1.  (x,z,X)  are  neglected  relative  to  Z.  U = t'Z 

3 9 

is  the  constant  sliding  speed.  Also  assumed  to  vanish  are  (ii2.V'?,t,  “r-,  ttt)  ■ 

3 2 f)  t 

The  differential  equation  is  reduced  to 


d ,Uph  1 , 3 ^ „ 

— ( o - — P ^)  = 0 


dx  ' 2 12p 


(4.1) 


The  boundary  conditions  are 


p (x  = 0,a)  = p^ 


(4.2) 


Non-dimensionalization  is  achieved  by  defining 


x = x/a;  h = h/c;  p = p/p„;  A = 


P C 
' a 


where  C is  a representative  gap.  A is  the  compressibility  number  for  a slider. 


Equations  4.1  and  4.2  are  then  rewritten  as 
— (Aph  - h^  p ^)  = o 


dx 


(4.3) 


p (x  = 0,  1)  = 1 


(4.4) 
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Equatlon  (4.3)  can  be  integrated  once  to  yield 


p-h  (1  - . h * 

dx 


— * 

h is  an  integration  constant;  it  represents  the  equivalent  gap  if 

profile  were  purely  of  the  couette  type.  Suppose  ^ vanislies  at  x 

dx 

* 

ii  = (pii)  I -* 


At  (x  = 0,  1),  p = 1,  then 


^ ^ (h  - I ) 

dx  ii 


Thus,  if  A is  to  be  recognized  as  positive  by  convention,  then  at 
edge,  the  sign  of  the  pressure  gradient  is  determined  only  by  (h  - 
following  observations  can  be  made. 


4.1.1  h converges  from  the  leading  edge.  An  extremum  of  [i 
exist  at  x where 


h = h (x  ) < ii  (0) 


This 
— * 

h < 


condition 


h (0);  p 


can  only  be 
>1 


satisfied  by 


Thus  pressure  near  a converging  leading  edge  is  ele\';i 
ambient  and  must  satisfy  the  ctnstraint  on  its  maxima 


1 < p*  = h 

b(x*) 


< b(0) 
il(x*) 


4.1.2  h diverges  from  the  leading  edge.  Using  similar  re  is> 
above,  one  concludes  that  pressure  near  .a  diverging  1 1, 
is  reduced  from  the  ambient  and  its  minima  is  hiaiki'ti 


/ ‘ '•nr'  ' '-*4  - 1 


*•  ^ 


(4.5) 

the  velocity 
*,  then 

(4.6) 

(4.7) 

either  ambient 
h'^).  The 

J. 

- p would 


ted  from  the 

(4.8) 

>n  i ng  as 
■ading  edge 
'd  in  t lie  range 

( t 


- I^^- 


h (0)  ^ 

h(x*) 


( '..9) 


The  above  results  ire  aflected  by  compressibility  mainly  in  the  lonstraints 
ot  pressure  extrema  according  to  (4.8)  and  (4.9).  Further  progr.  as  va)iild 
depend  on  specific  details  of  h.  Of  special  interest  is  the  step  slider, 
shown  in  Figure  4.2,  which  can  be  solved  in  close  form  ll2). 


4.2  Infinitely  Long  Step  Sliders  (Figure  4.2) 

h is  sectional Iv  constant  in  the  two  regions  (0  ' x ■ x , x < x I). 

o o 

Equation  (4.5)  can  be  rewritten  as 

^ (I  - ^)  (4.10) 

dx  p 

where  y = h /h  and  ‘ = //h  are  constants  in  eacli  region. 

From  Equation  (4.10),  one  re.'.sons  that  p must  be  monotonic  in  the  range 
0 X X ■'  x^.  If  one  further  imposes  Equation  (4.5)  as  ;i  matching  condition 
at  the  step, 


h* 


h {1  - (x^  )} 

dx  ^ p(x  ) 


1 - (X  +) 

dx  " 


then  one  would  reach  the  conclusions,  for  h^  ' 1 


1 < 


P(x^) 


' h ^ ^ 


(9.11) 


(4.12) 


in  order  to  satisfy  the  boundary  conditions  at  (x  = 0,1).  Equation  t.4.10) 
can  be  formally  integrated  for  either  region  to  yield; 

P-^ 

p - 1 + Hn  (-j— T-)  = )ijX  (4.13) 


for  0 < X < x^,  where  = h*/ii^  and  = A/h^;  and 


p - 1 + h*  nn  (^ — 2.)  ^ A(x-l)  (4.14) 

h*-l 


f or  X < X < 1 . 
o 


These  equations  must  result  in  the  same  p(x  ) = 


— — — * 

P h -h 

h (p  -1)  + h*  ?n  (— °—  ■ ■ — ) = — X 

^ ^ h,-h*  h^ 


p^,  llaretore  one  obtains 


a . I b ) 


h*-p 

p -1  + h*  Hn  (— -)  = A(x  -1) 

o h*-l 


(4. lb) 


They  define  simultaneously  the  numerical  constants  (j'^,  h ) for  eicli  set  of 
(hj , x^;  A).  The  correct  solution  can  be  found  in  the  raiiRP 


1 < p < h*  < li, 
*^o  1 


(4.17) 


At  extreme  values  of  A,  the  limiting  behaviors  can  be  examined  analytica 

4.2.1  Small  A (Incompressible  Case) 

One  mav  seek  the  solution  in  the  f(?rm 


1 + Att  + 01  A") 


(4.18) 


Substituting  into  Kquations  (4.15)  and  (4.16),  collecting  terms 
in  the  first  power  of  A,  <ine  finds  two  simple  linear  algebraic 
equations  which  can  be  readily  solved  to  yield 


r (1-x  )h,+x  ^ 

I r-—  h 71 

/ 1 “ \ r 3 . “ 3 

. ( l-x  )h,+x 
o 1 o 


x (l-x  ) (li  -1  ) 
o q 1 

(l-x  )h^  + X 
o 1 o 


(4.19) 


Accordingly, 


p(0<x-  x)  = I+Ati 

- o o - 

X 


X 


o 


p(x  <X<1)  = 1+Ati  (---—) 

o - - o , ■ 

1 -X 

o 


'’LL 


f - ' '-*d  1 


(4.20) 
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These  seem  to  be  the  incompressible  step-b(*nr  i iik  solutions. 

4.2.2  Lnrge  A (iliglily  Compressibie  C.-ise) 

i'he  limiting  behaviors  of  Equations  (4.16)  and  (4.16)  for  large 
h require 

Jim  (P|^>  ^ ~ (hj  , hj ) 

Writing  6h  = hj  - h*  and  <5p  = h^  - p^;  tlien  Equations  (4.15)  and 
(4.16)  become: 

' h . (h  -1  )+'Sh-h  j {p 

h (h  -l-5p)  + (h,-4h)  /n  7 “ (4.21) 

1 ^ ‘ 6h  h ^ 


h,-l-4p  + (h  -6h)  dn  — ^ — 

^ ^ (hj-l)-6h 


•(1-x  ) 
o 


For  large  A,  dominating  terms  for  the  respective  lefthand  sides 
of  the  above  equations  arc  those  containing  tn6h  and  C,n(6p-eh). 
Thus,  for  a first  approximation. 


Ah  ' _ h.E..  (1-E,,) 

- I 

1 + E ' 

" I 

Ap  t ,h  E (1-E  )+E  (1-E  ) , 


(4.23) 


where  (E  , E,  , E~)  are  exponentially  vanishing  functi('ns: 
o 1 Z 

I Ax^. 

E,  = exp  I rj  ; E = exp 

1 . I., 

E = E-{(h,-l)  - h.E  } (4.24) 

o 11  1 Z 

Accordingly,  upon  substituting  Ah  into  Equation  (4.13)  and 
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retaining  only  terms  with  commeasurate  accuracy,  one  finds 


p(0<x^x^^)  = 1 -t- 


1+E 


' exp  - 


H,' 


p(x^  < X < 1)  = 1 + 


(h^-D(l-E^) 

r+E 


1-exp  - 


A(]-x_)^ 
h ' 


(4.25) 


In  essence,  p has  become  almost  sectional ly  uniform; 

p(0  < X < X ) = 1 and  p(x  < x < 1)  = h,,  except  for  boundary 

t>~_  i 

layer  behaviors  as  x ‘ (x^  , 1 ) ■ 

Sample  results  for  (hj^  = 1.5,  x^  = 0.  50)  are  shown  for  a full 
range  of  A in  Figure  4.3.  The  primary  effects  ot  .ompressi- 
bility  are  related  to  the  pressure  peak  (which  occurs  at  the 
step  junction)  and  the  pressure  profiles  on  its  two  sides. 

For  small  values  of  A,  the  amount  of  peak  pressuri'  elevation 
above  the  ambient  is  essentially  proportional  to  the  product 
of  the  sliding  speed  and  the  fluid  viscosity  and  is  insensitive 
to  the  ambient  pressure  level.  At  large  compressibility 
numbers,  the  peak  pressure  asymptotically  approaches  a constant 
multiple  of  the  ambient  pressure;  this  asymptotic  multiple  is 
numerically  equal  to  the  gap  ratio  h^ . The  pressure  profile  is 
approximately  triangular  at  a small  compressibility  number, 
which  causes  concavity  in  the  inlet  portion  and  convexity  in  the 
exit  portion  of  the  pressure  profile.  As  A becomes  asvTnptoti- 
cally  large,  the  pressure  profile  approaches  two  horizontal  lines 
with  a rapid  increase  at  the  step  and  a rapid  decrease  at  the 
exit  edge. 


4.3  Infinitely  Long  Taper  Slider  (Figure  4 ._4_) 

For  a taper  slider,  the  film  thickness  may  be  described  as 

h = 1 + (1  - x)  Ah  (4.26) 


- J9- 


so  that 


dh 

dx 


- Ah 


(4.27) 


which  Is  a constant.  This  fact  allows  one  to  reduce  Kquation  (4.5)  into  an 
exact  integral  with  the  aid  of  the  transformation 

Ip  = ph;  df  = - ~ dx;  A = P.nh  (^.28) 

h 

Equation  (4.5)  accordingly  can  be  rewritten  as 

T = h*;  > = 2!L  (4.29) 

Ah 


This  equation  can  be  integrated  in  closed  form;  whereupon  substitution  of 
the  original  variable  h for  inf  yields 

X 

(^-s^Xu^-e^)  1-B2  ■i'-s, 

(l-Bj^)(l-B2)  ^ ^1^^  ^ 30) 


where 


(81,  Bp  - 

X ^ ' X , 

2±  ; T - 

(4.31) 

For  small  values 

of  A,  (e^,  $2)  are 

comp  lex : 

(6^.  B2)  = 

1 ± iY;  Y = V 

>h*-A 

(4.32) 

Then,  instead  of 

Equation  (4.30) 

^ (<i-  - |)^  + 

2L1  'a.  -i_ 

y(’p-l)  ' rZ 

(4.33) 

« - . 

ty  6xp  Y wdn 

. (1 
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For  each  X , h is  found  by  satisfying  the  inlet  boundary  condition  (i|/  = i + Ah 
when  ii  = 1 + Ah)  in  citlier  Equation  (4.30)  or  (4.33),  de[)ending  on  whetiicr 
(Ej^,  62)  ate  real  or  complex.  This  complicates  the  computation  process  some- 
what; the  general  numerical  features  are  as  sliown  in  Figure  4.5.  At  small 
A,  p is  directly  proportional  to  A so  that  the  actu/il  lilm  pressure  is 
independent  of  the  ambient  pressure.  As  A becomes  increasingly  larger, 
however,  p does  increase  indefinitely,  but  becomes  ultimately  bounded  by  the 
asymptotic  limit  ph  = 1 + Ah  whereupon  at  tlie  trailing  i^dge  p drops  abruptly 
to  meet  the  ambient  condition. 


4.4  Plain  dournal  Bearing  (Figure  4.6) 

The  simplest  journal  bearing  is  a rotating  journal  with  its  axis  displaced 
from,  but  remaining  parallel  to  the  axis  of  the  stationary  bushing.  Both 
the  journal  and  the  bushing  surfaces  are  purely  cylindrii  al.  Axial  motion 
is  assumed  to  be  absent.  The  following  conditions  thus  prevail  '‘or 
Equations  (3.9)  and  (3.10) 

= w;  ^2  = 0 


h = C(1  - CCOS0)  (4.34) 


Using  the  dimensionless  variables 


z = z/R;  P = P/Pg'  h - h/C  = 1 - ccosO 


A 


6uuiR 
Pa  ^ 


2 

2 


(4.35) 


then  Equation  (2.24),  for  impermeable  walls  and  at  steady  state,  is  reduced 
to 


9 ,r3-  2P\  . 
— (h  p -r^)  + 


90 


90' 


9z 


(h  p 

9z 


(ph) 


(-).36) 
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The  boundary  conditions  are 

= + L/D)  = 1 (4.37) 

and  periodicity  in  0. 

In  its  complete  form.  Equation  (4.36)  is  not  amenable  to  a rij;orous  clo.sed 
form  solution.  Results  obtained  by  numerical  computation  have  been  quite 
thoroughly  compiled  for  engineering  analysis  fl3,14].  Asymptotic  approxi- 
mations have  in  some  instances  rendered  closed  form  solutions,  while  in 
other  cases  simplified  the  numerical  formulation  considerably.  Ihese  are 
reviewed  in  the  following  to  bring  out  the  relevant  plivsical  features  which 
are  substantially  valid  irrespective  of  the  asymptotic,  approximations. 

4.4.1  Small  Eccentricity  Solution 

Presuming  analytical  continuity  in  the  limiting  condition  ' ' 0, 

a power  series  solution  may  be  postulated: 

p = 1 + cp  (4.38) 

E 

Substituting  into  Equation  (4.36)  and  collecting  term.s  associated 
with  the  first  power  of  r,  one  obtains 

. 2 

^ l(^  - A)p^}  + -^  = /.sin0  (4.38) 

3 z 

and  the  corresponding  boundary  conditions  are 

p^(z  = + L/D)  = 0 (4.40) 

This  formulation  is  knowm  as  the  "l.inearized  P"  solution  [13,161 

2 

and  is  valid  for  all  (A,  L/D)  provided  o{ < ‘ is  sufficiently 

small . 

An  analogous  formulation  can  be  carried  out  with  the  expansion 


il)  = ph=l  + fil^  (4.41) 
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Accordlngly,  Equation  (4.36)  is  reduced  to 


Je  ^ ■ 2 


= COS0 


(4.42) 


with  the  boundary  conditions 


<P^(z  - + L/D)  = - COS0 


(4.43) 


This  latter  formulation  is  known  as  tlie  "Linearized  PH" 
solution  [17]. 

It  is  almost  self  evident  that  n and  i are  related  to  each 

E 

Other  by 


sj;  = p - cos6 
e e 


(4.44) 


The  advantage  of  treating  <p  instead  of  ph  as  tlie  unknown  is  duo 
to  the  asymptotic  condition 


llm  ip 

\-KX, 


= constant 


(4.43) 


for  all  E,  even  though  the  constant  itself  is  si'mewhtt  depemient 
on  e.  This  interesting  result  v.;i11  be  further  elaborated  in 
Section  4.4.2.  For  any  finite  A,  however,  does  not  have  any 
advantage  over  p^.  In  the  present  section,  no  further  attention 
will  be  given  to  ip  . 


The  harmonic  dependence  on  0 in  the  nonhomogeneous  term  of  Equa- 
tion (4.39)  and  the  inherent  0-periodicity  of  p^  suggest  the 
solution  to  be  of  the  form 


Re{ qe ' ^ } 


(4.46) 


Accordingly,  Equation  (4.39)  is  reduced  to 


{— 7T  - (1  + iA)}q  = - iA 


(4.47) 


IT-  ' '-*4  - 
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A closed  form  solution  satisfying  tlie  boundary  cond  tion 

q(i  = + L/D)  = 0 (4.48) 


^ - ^iTIa>  ■ ’ 


cosh  v'  1+iA  z 


cosh 


. (^)J 


(4.49) 


Since  p as  given  by  Equation  (4.46)  is  a sinusoidal  function  of 
0,  real  and  imaginary  portions  of  q represent  respectively  a 
component  in  phase  with  the  minimum  gap  region  and  a quadrat ive 
component  shifted  a quadrant  toward  the  convergent  side  of  the 
bearing  gap. 


, . . , r 1 coshz  I 

lim  q = lAl  1 = — } 

,A->0  cosh(p) 


(4.50) 


1 im  q = 1 - 


cosii 


. ii\  z 


cosh  iA 


(--) 


(4.51) 


Thus,  at  small  A,  q is  proportional  to  A and  is  a purely 
quadrature  component.  At  large  A,  except  for  the  end  regions 
:iA  (L/D- I z 1 ) = 0{ 1 } , q is  essentially  uniform  with  respect  to 
z and  is  inphase  with  the  minimum  gap. 

4.4.2  Large  A Solution 

For  large  A,  Equation  (4.36)  may  be  rewritten  as 


30 


1 r 3 ,r3-  . 3 ,r3-  , 

3z  3z 


(4.52) 


where  ip  = ph.  If  one  presumes  the  terms  enclosed  in  the  brackets 
to  be  bounded,  then  as  A -*■  Equation  (4.52)  would  indicate 


lim 

A->“> 


3(ij) 

30 


= 0 


(4.53) 


s . * -■ 

• /•_  .1  "• 
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This  condition  is  a bit  unusual  in  two  aspects.  Firstly,  Equa- 
tion (4.53)  is  in  direct  contradiction  witli  the  boundary 
condition 


86  - 

■ z=+L/D 


Secondly,  Equation  (4.1)5)  has  a total Iv  different  structure  as 
compared  with  Equation  (4.36)  from  the  standpoint  of  the  tlieory 
of  partial  differential  equations. 

To  remove  the  above  two  objections,  one  must  reformulate  tlie 
present  problem  in  terras  of  a "matched  asymptotic  expansion" 

[18].  With  such  an  approacli.  Equation  (4.53)  is  to  In-  imposed 
only  in  the  open  domain 

|z|  ^ ^ (4.55) 


as  an  "outer  solution."  Whereas,  in  order  tliat  tlu>  boundary 
conditions  may  be  reconciled,  an  edge  region  bordering  each 
end  would  be  recognized  such  that 


= Of,  A'  } for  0 < - |z!  = 0{"  -'1 

8z 


(4.56) 


Accordingly,  a rescaling  of  the  axial  coordinate  is  in  order; 
namely,  for  the  edge  regions,  one  can  define 


A'  (jt  - 1^1) 


(4. 57) 


so  that  instead  of  Equation  (4.52),  one  lias  the  edge  equation 


^.2^2^  ae 


where  Neglecting  the  righthand  side  of  Equation 
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(4.58)  no  longer  causes  ttie  two  objections  previously  raised 
with  respect  to  Equation  (4.53).  The  boundary  condition  at 
4=0  can  be  directly  identified  with  Equation  (4.54). 

The  procedure  to  connect  the  "outer  solution"  witli  tiie  boundary 
layer  solution  remains  to  be  developed.  The  key  is  contained  in 
the  "mass  content  rule"  for  self-acting  gas  bearings  [4]  wliich 
is  to  be  derived  in  the  following. 

Because  ii  does  not  depend  on  z,  Equation  (4.36)  ran  I'e  rewritten 
as 

“ f.' (ph)  - h^p  (4.59) 

3z 

Integrating  both  sides  over  the  full  <•  i reumt erence , as  p and  h 
are  both  periodic,  the  righthand  side  would  vanish;  so  that 

,2  '2-n 

(h  p )d0  = 0 (4.60) 

dz  - 0 


or. 


'2ti 

0. 


(h^p^)de 


(4.61) 


Since  p(z  = + L/D)  = 1,  one  obtains 


= 0 


-2ti 


'0  j, 


ii^de 


2ti(1  + -I  “) 


Therefore 


2n  Jo 


-3-2 

(h  p )dfl 


(4.62) 


This  is  the  celebrated  "Mass  Content  Rule. 
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Equation  (4.62)  is  valid  for  all  z.  If  one  applies  it  to  tlie 
outer  solution,  then 

2 1 ' r m 

lira  lidO 

2 3 2 

= limij/  = l+-:7f 


Or, 


lim  'p 

A'VCO 


I im  p 
A-x*' 


2 


e 


1 - CCOS0 


('i.h3) 


As  the  "edge  solution"  must  blend  into  the  "outer  solution" 
smoothly,  this  also  becomes  the  second  boundary  condition  for 
Equation  (4. 38) : 


P , (4- 

edge 


,e) 


3 2 ■' 

1 - > cosQ 


(4.64) 


It  is  clear  now  that  ; = ph  • constant  is  a condition  prevailing 
at  large  A for  all  ..  The  asymptotic  value  of  however,  is  an 
explicit  function  of  ‘ according  la  Equation  (4.63).  The  global 
inaccuracy  of  this  approximation  is  Ofl/i/A  which  correspomis 
to  the  axial  fraction  of  the  edge  region.  Thus  the  "linearized 
PH"  solution  actually  does  not  accomplish  any  improvement 
beyond  = constant  in  the  central  region.  The  term  given  by 
Equation  (4.44)  furnishes  a portion  of  the  next  improvement, 
which  is  of  0{1/,'T'}  in  the  global  sense. 


To  extend  the  accuracy  in  the  large  A region,  the  outer  solution 
may  be  found  by  solving  Equation  (4.52)  with  the  expansion 


r 
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3 2 - 1 

i()  e + A iji  -t-... 


(4.65) 


obtaining 

4 ■ 


*10  - 1 * ‘ * i ' ' 


is  an  additional  constant  to  be  determined  bv  Equation 
(4.62)  and  turns  out  to  be 


'^10  = « 


Thus , 


e 1 + -17  * ^ 


(4.66) 


The  "edge  solution"  should  be  similarly  exi>anded  as 


it  - -I-  A ij  . + 

^eO  el 


(4.67) 


The 


leading  terms  are  to  satisfv  the  followinj'  conditions. 


3ij) 


3G 


- h dr  (* 


3i(i 


eO 


34  '^eO  34 


) = 0 


(4.68) 
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ij^^^(c=o,e) 


3 3 2 

I r 1 + - E .sine 


Since  further  correction  would  be  of  Of  A for  bt)t  li  the 
"outer"  and  the  "edge"  solutions,  loi al  accuracy  up  to 
0{A  is  achieved  by  (ifiQ,  ever;/where.  rlobal 

accuracy,  however,  is  of  0{A 


4. A.  3 Long  .Journal  Bearing 

If  L/D  is  very  large,  (L/2)  instead  of  R is  a inorc  .appropri- 
ate scale  for  the  axial  coordinate.  If  one  defines 


n = 2z/L  = (7)  z 


then  Equation  (4.36)  becomes 


^ /r3-  ^ 

( }■»  n ^ — An 


Aph)  = 


,D.2  1 ,p3-  ep, 

- 3^  P -.V 


Clearly,  if  one  lets  the  righthand  side  vanish,  one  again 
obtains  an  "outer  solution"  which  is  not  valid  for  I - |n|  = 
0{D/L}.  The  outer  solution  is  governed  by 


^ ^ - Aph)  = 0 


Or,  integrating  once. 


-3-  — * 

h p dQ  = Uph  - 'I'  ) 


)|;  is  a constant  of  Integration.  The  value  p^  = p(fl=0)  in.av  he 
regarded  as  the  second  constant  of  integration,  tinis 


p(0)  = p + A d0'  (3^  - ^^^) 

0 h h p 


/ ' 


V • - -V 


r 


For  every  i|i  , the  correct  value  of  is  found  l)y  Llie  condition 


p(+2mr)  = p(0) 


(4.74) 


A suitable  continuous  iteration  procedure  is  to  increment  n 
repeatedly  to  find  p(0): 


Pc  = Pq'  ^ ^ 


i 1 ih  , 

-0  h h p 


(4.75) 


Integration  in  the  negative  (J  direitiun  assures  r.umi  r i ca  1 
stability.  Convergence  is  reached  when  furttier  incri'iiient  of 
n cease  to  yield  significant  change.  WTieriupun  t ne  "mass 
content  rule"  would  be  checked  to  verify  the  lorrectness  of  , 

The  ambient  boundary  condition  is  satisfice  by  the  v)nstruc- 
tion  of  an  edge  solution.  The  a.\ial  coordin.itc  for  tlic  edge 
solution  is  to  be  rescaled  in  order  to  ret.ain  the  righthand 
side  of  Equation  (4.71).  Thu.s  tiefjning 


(^)  (1 


l|) 


(4.76) 


at  either  end,  the  edge  problem  is  expressed  as 


^ (h3  - I"-- 

30  Pe  39 


J -3  - ^Pe 

A Pe^)  H (h  p^  -J)  = 0 


(4.77) 


The  required  boundary  conditions  include 


P^(9+2ti,c)  = p^(0,C); 


30 


3p„ 

= —^1 
-3^  1 

0+-2ti,5  |0,f. 


p^(0,r,=O)  = 1 


p^(0,C-^)  = p(0) 


(4.78) 


L 
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The  last  expressions  ensure  smooth  blending  of  into  the 
"outer  solvition." 


4.4.4  Short  Journal  Bearing 

If  L/D  is  small,  the  compressibility  number  would  have  to  be 
rescaled  as 


A 


sh 


A 


Accordingly,  Equation  (4.71)  can  be  rewritten  as 


(4.79) 


d 


<■^'5  - 


-JL 

'sh  30 


(ph) 


(4.80) 


The  righthand  side  is  dropped  when  (—)  <■  0.  If  n is  replaced 

by  > Equation  (4.80)  becomes  identical  as  Equation 

(4.58),  which  is  concerned  with  the  edge  problem  with  a large 
compressibility  number.  The  axial  boundary  conditions  are, 
however,  different.  Here,  one  has 


p(n  = + 1 , 0) 


(4.81) 


.■ 


4.5  Nomenclature 


Symbol 

a 


c’o.(-i 


E E K 
0’  *^1’  2 


h 

h 


_* 

h 


4h 


Ah 

L 

P 

Pa 

Pe 

P 

Pe 

Po 

_* 

P 

6p 


q 

R 

u 


Def Inlt ion 

Slider  width  in  x-direction 

Constants 

Journal  diameter 

Rxpoential  functions  in  compressible  step  slider 
solution,  eq.  (4.24) 

hear  in;’,  nap,  h,,  - hj 

h/C 

Inlet  to  exit  gap  ratio 

C('iist  int  of  intenration  in  slider  solution 

- h for  step  slider  at  large  A 

h I - 1 of  taper  slider 
Journal  length 
Fluid  pressure 
/\mbi(>nt  pressure 

Perturbation  of  p with  respect  to  t 
p/p., 

Edge  solution  of  p for  long  journal 

p at  X of  step  slider 
o 

Extremum  of  p 

iij  - p^^  of  step  bearing  for  large  A 

Function  in  solution  of  p^  for  Journal,  eq.  (4.46) 
Journal  radius 
Sliding  speed 


Chapter/ 
Sec  t ion 


4.2 


4.2 

4.3 


' V*- 


• ^ 
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Svmbol 

X 

X 

X 

o 


Pet  Ip  it  ion 

Coord ina  Le 
x/a 

Width  fraction  of  step 


(Mi.ipLer/ 
Sec  t i on 


_* 

X 

z 

z 


Y 

Y 


e 

4 

n 

e 


0’ 


\ 

> 

(li 

A 

A 

A 

sh 


;i 


Value  of  X at  pressure  peak 
Axial  coordinate  of  iourna] 
z/R  of  journal 

Numerical  parameters  in  taper  slider  solution, 
eq.  (4.31)  and  (4.32) 

—•k  — 

h /h  for  step  slider 

Numerical  parameter  in  taper-slider  solution, 
eq.  (4.32) 

Eccentricity  ratio  of  journal,  e/C 

Dimensionless  axial  edge  coordinate  for  long  journal 

Dimensionless  axial  coordinate  of  long  journal, 

2 z / 1, 


Polar  angular  coordinate 

Dummy  variable  of  Integration  for  0 

A/h'^  for  step  slider 

A/Ah  for  taper  slider 

Rotational  speed 

Compressibility  number  for  slider. 


-^-'4 

Pa^' 


Compressibility  number  for  journal. 


6pif  R‘ 
Pa^"  ' 


Short  journal  compressibility  number,  (— ) A 


Viscosity  coefficient 

Independent  variable  of  taper  sJider  solution  fnh 


4.3 

4.2 

4.3 

4.4 
4.4.3 
4.4.3 


4.2 

4.3 


4.1 

4.4 

4.4 


4.  3 
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Svmbol 

n 

o 

Ip 

'Pi 


Ip  1 Ip  , 
’^eo  el 


Definition 

Incompressible  step  slider  solution 

pii 

First  order  term  in  reciprocal  expansion  for  the 
"outer  solution"  of  journal  at  larpe  A,  eq.  (4.66) 

"Edge  solution"  of  ip  for  journal  at  large  A 

Leading  terms  in  reciprocal  expansion  for  ip 
eq.  (4.67) 

Perturbation  of  ip  with  respect  to  c 

Constant  of  Integration  in  solution  for  an 
infinitely  long  journal 


Chapter 
Sect  ion 


4.2 


4.4 


4.4 


4.4 
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5.0  EXTERNALLY-PRESSURIZED  BKAKINC .;S 

In  an  externally-pressurized  gas  bearing  the  pressurization  level  (relative 
to  the  ambient)  set  the  range  of  film  pressure.  Due  to  compressibility,  the 
supply-to-ambient  pressure  ratio,  instead  of  the  pressure  difference,  is  the 
governing  parameter  in  the  operating  characteristics  of  the  bearing.  The 
steady-state  governing  equation  for  the  most  general  externa  1 1 y-pi essur i zed 
bearing  is 

-*■  1 >■  . 1 2u 

'2  ' (h  p 2 p)  = 7—  Hi  (5.1) 

f accounts  for  effluxes  through  both  bearing  surfaces.  For  exieniallv 
pressurized  porous  gas  bearings,  >' ! wtnild  normally  be  negative. 

In  comparison  to  the  self-acting  gas  bearing,  mathematics  reqiiireJ  to  st'lve 
the  externally  pressurized  gas  bearing  is  relatively  straight  forward.  The 
unknown  of  Equation  (5.1)  may  be  considered  to  be  ~ and  the  enuation  is 
in  fact  linear.  In  many  instances,  whet,  the  gap  is  unitorm  and  tlie  surlaci-s 
are  impermeable,  closed  form  exact  solutions  are  avail.ible.  Hov.ever,  to 
analyze  an  externally-pressurized  bearing,  most  of  the  time,  one  cannot 
treat  Equation  (5.1)  as  .in  isolated  problem  because  Liie  bearing  film  is 
often  only  a portion  of  the  overall  pneumatic  network  between  the  pressur- 
ized source  and  the  ambient.  This  is  because  an  externa  1 1 y-ju  es^u r i zed 
bearing  generally  wmild  require  some  kind  of  upstream  restriction  to 
realize  a stiffness.  it  is  when  the  restriction  system  is  incliuU'ci  that 
the  overall  analysis  can  become  substantially  more  complicated. 

in  a porous  bearing,  restriction  takes  place  in  a distributed  manner.  Its 
complete  analysis  involves  solving  the  pressure-flow  field  problem  Logetiier 
with  Equation  (5.1)  with  appropriate  pressure  and  flow  continuity  conditions 
at  the  interface  between  the  porous  material  and  the  bearing  film.  The 
porous  bearing  surface  itself  introduces  some  additional  complication  since 
the  non-slip  condition  on  the  velocity  profile  would  have  to  he  revised. 

Some  of  these  questions  are  topics  of  current  research. 

For  bearings  with  impermeable  walls,  there  are  three  major  discrete 


restrictor  types;  namely,  the  capillary  restrictor,  the  orifice  restrictor, 
and  the  inherent  restrictor.  The  most  commonly  used  restrictor  type  is 
the  last  one — the  capillary  restriction  seldomly  can  be  achieved  with  a 
practical  configuration,  the  orifice  type  restrictor  is  vulnerable  to 
"pneumatic  hammer”  instability. 

In  the  following  sections,  general  features  of  the  inlierent  restrictor  will 
first  be  discussed;  tlten  representative  solutions  of  Equation  (5.1)  for 
impermeable  walls  will  be  presented  with  allowance  for  the  use  of  inherent 
upstream  restriction. 

5.1  The  Inherent  Restrictor 

When  gas  enters  into  the  bearing  gap  from  an  essentially  stagnant  state  in 
the  supply  reservoir,  its  pressure  undergoes  an  abrupt  drop  due  to  the 
formation  of  the  velocity  head.  Subsequently,  under  the  action  of  wall 
shear  stresses,  the  velocity  profile  gradually  develops  into  a parabolic 
distribution  whereby  its  pressure  undergoes  an  additional  drop.  The  flow 
length  required  for  the  velocity  profile  development  process  to  complete 
is  proportional  to  the  flow  Reynolds  number.  Due  to  these  phenomena,  the 
effective  film  inlet  pressure  is  lower  than  the  supply  stagnation  pressure 
by  an  amount  which  is  proportional  to  the  velocity  head  at  film  inlet. 

Thus,  if  the  supply  hole  is  properly  sized,  to  make  this  pressure  drop 
significant,  there  is  an  "Inherent  inlet  restriction." 

A theoretical  model  for  the  behavior  of  the  inherent  restrictor  is  the 
entrance  of  an  incompressible  viscous  flow  into  a parallel  channel  [19,20j- 
According  to  this  calculation,  the  static  pressure  distribution  at  the 
channel  entrance  can  be  described  in  a universal  manner  as  shown  in 
Figure  5.1.  The  symbols  used  in  labeling  the  graph  are: 

p^  static  pressure  at  entry  x distance  from  entry 

p local  static  pressure  v kinematic  viscosity  of  liquid 

p density  of  liquid  h channel  height 

U mean  velocity 
m 
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The  static  pressure  distribution  has  an  asymptotic  tangent  which  is  in 
agreement  with  the  lubrication  theory  (fully  developed  viscous  channel 
flow).  According  to  this  result,  if  one  replaces  the  curve  in  Figure 
5.1  by  its  asymptotic  tangent,  the  effective  entry  film  pressure  is  Iowit 
than  the  upstream  stagnation  pressure  by  the  amount 


pu; 

Po  - Peff  = 0-^0^  > 


(5.2) 


If  one  wishes  to  account  for  the  deviation  between  tlie  asymptot  io  L.ingeni 
and  the  actual  curve,  a "displacement"  lengtii  can  be  calcul.ited  lo  lie 


6 


+ 


0.00623 


(5.  ij 


The  integrated  effect  of  this  displacement  length  is  a load  del  ic  iimcv  of 
the  amount 


6W  = 4.692  X 10 


1 


(5.4) 


This  is  a second  order  effect  which  is  much  smaller 
neglected  in  the  lubrication  theory.  Therefore,  it 
Kq.  (5.2)  to  aircount  for  a restrictor  pressure  drop 
regarding  the  details  of  pressure  development  in  an 
film. 


than  the  terms  which  an 
is  sat isf ac I ory  to  use 
without  further  cinii  ern 
entrance  rcgi<'n  ot  the  t luid 


The  asymptotic  law  for  inherent  inlet  pressure  drop  ;ts  given  by  Fq.f5.21  is 
mainly  accountable  by  the  nondissipat ive  conversion  of  pressure  to  kinetic 
energy  [20].  Therefore,  one  may  extend  this  idea  to  compress  i li ! e fluid  dvnaniiis 
as  a basis  for  calculating  the  pressure  drop  through  ;ni  inherent  resii  ictor. 

The  following  conceptual  steps  are  assumed  in  order  to  develop  the  applic.tble 
computational  procedure  (see  Figure  5.2): 

• The  bearing  wall  is  in  thermal  equilibrium  with  tlie  reservoir  gas. 

■ A thcoat  section  is  recognized.  Since  the  inlet  edge  in  a practical 
construction  would  be  more  or  less  rounded,  the  thro.u  st'ctii-cn  mav 
be  defined  as  the  termination  point  of  the  edge  round. 


-(.J- 


The  flow  process  from  the  reservoir  to  the  throat  follow  the 
iseiitropic  law. 

Near  the  throat  region,  heat  transfer  between  the  walls  and 
the  gas  flow  takes  place  rapidly  to  equalize  their  tempera- 
tures with  little  pressure  change. 

Further  development  of  the  velocity  profile  involves 
negligible  compressibility  effects. 


Subscripts  "s",  "th",  and  "in"  will  be  used  to  designate  the  respective 
states  of  supply  stagnation,  throat  and  effective  film  inlet. 


Let  the  mass  flux  per  unit  inlet  edge  l^e 


^ = ‘^th^hh 


According  to  the  isentropic  law. 


"s  Ps 


T V-I 

th  _ Pth.  7 

T ' “ p ^ 

s s 


0-1) 


Y-1  ^ 2 


(3.6) 


Thus,  Equation  (5.5)  can  he  rewritten  to  express  the  mass  flux  in  terms  of 
the  stagnation  states  of  the  reservoir  and  the  pressure  ratio  (P(.},/Pj.O 


Y-1-1/2 


2Y  ,^th.- 

(Y-1>RT  Ps  p ^ 
s s 


th,  Y 


With  temperature  equalization  (the  pressure  being  unchanged),  the  inlet 
velocity  head  becomes 
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i ( -t/h) 

s 


5^T  <r 
± s 

2 u2 

Pth^ 


Now,  making  use  of  Equation  (5.2) 


p . = p.  + 0.601  pU  ) . 

th  '^in  2 m in 


C5.9) 


Combining  Equations  (5.7),  (5.8),  and  (5.9),  one  obtains 


Pin  Pth  7 ■ ^s  Y ^th  Y 

_\n  = ^ _ 0.601  (~)  i (-^)  - (--y  (-^'i)'^ 
P..  P.-  P,i,  P,a.  P., 


th  ' th 


(5.10) 


For  gas  lubrication  calculations,  it  is  more  convenient  to  use  the  pressure 
flux  \l>  instead  of  the  mass  flow  <1>;  that  is 


>p  = " T $ 


so  that  Equation  (5.7)  can  be  rewritten  as 


p . ‘ T h 
s s 


Y-1 

(^)  " 

Ps 


(5.1!) 


Equations  (5.10)  and  (5.11)  allow  one  to  calculate  p.^  upon  given  the 
stagnation  state  variables  of  the  reservoir  and  the  pressure  flux. 


Since  the  function  of  the  restrictor  is  mainly  in  the  preservation  of  an 
equilibrium  condition,  it  is  essential  to  know  the  differential  perturba- 
tion characteristics  of  p^^  with  respect  to  p and  h.  This  is  obtained  !\v 
differentiating  Equations  (5.11)  and  eliminating  the  differentia!  of  Pfi/Pj.- 
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Craphical  presentations  of  Equations  (5.10),  (5.11),  and  (5.12)  are  shown 

in  Figure  5.3.  Despite  these  rather  complicated  expressions, 
the  following  observations  can  be  noted: 

• Compressibility  affects  the  restrictor  pressure  drop  only  modestly. 

Even  if  the  flow  should  be  sonic  at  the  throat,  the  effect  of  inlet 
flow  development  is  within  5%  of  the  incompressible  case.  It  is  higher 
for  a monatomic  gas  but  lower  by  a similar  fraction  for  a diatomic  gas. 

• The  relationship  between  the  dimensionless  flow  rate  and  the  normalized 
supply-to-throat  pressure  drop  is  approximately  an  elliptical  quadrant 
for  both  monatomic  and  diatomic  gases. 

• The  pressure-flow  perturbation  relationship,  when  expressed  in  a 
dimensionless  representation,  is  essentially  same  for  both  monatomic 
and  diatomic  gases. 

• Perturbation  of  restrictor  pressure  drop  by  flow  is  relatively  insig- 
nigicant  unless  the  supply-to-throat  pressure  is  at  least  50%  of  the 
critical  pressure  drop. 

Based  on  these  observations,  the  following  empirical  formulae  are  proposed. 


Pth  Pii 


0.601^ 


1 + C„( 


Ps-Ptli 

Ps~Pcr 


T 

) (Ps  - Pth) 


v =-A/(P§Jlth)(2  - £s:JPt]l) 


Ps  Per 


Ps  Per 


(5.13) 


P(,j.  is  the  critical  pressure  for  sonic  flow,  is  the  corresponding  flow 

function  obtained  by  substituting  p^-j-  for  pj.j^  in  Eq.(5.11)  and  is  an  empirical 
numerical  constant.  Their  numerical  values  are 
Monatomic  Diatomic 


p /p  0.48714  0.52828 

*^cr  s 

0.72618  0.68473 

Cq  0.1139  -0.1150 

These  empirical  formulae  give  "exact"  results  at  either  extremes  of  Pj,^-^p^jj<^Pg . 
At  the  midpoint,  Pj,j^  = 1/2  (Pj,j.  + Pg)  > the  inaccuracies  arc 
Error  Monatomic  Diatomic 


% (Pth-Pin)  1-2^^ 

% f 0.496 


0.742 


wiilch  are  remarkably  small. 


j*. 


-0.139 
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In  Eqs.(5.13),  Pg  - can  be  eliminated  to  render 


0.601V1  + 1 X 

^cr  '' 


where 

X = (is . f-tji) 

Ps  ' Per 


(5.14) 


Performing  a variation  operation  on  Eq.(5.14),  it  is  found 


6p 


in 


(Ps  ■ 


) ] 
; r 


1 + 0.601 


(1  + 2CoX2)'| 


+ CqX 


r 


^erV^ 


(5.15) 


Since  it  has  been  shown  in  Figure  5.3  that  restrictor  pressure  drop  is 
most  likely  to  be  at  least  half  of  that  for  sonic  flow,  one  may  safely 
approximate  C^X  by  in  these  formulae.  Thus  one  obtains 


Pin  Ps  ^Ps  ~ Per^ 


1 + 


O.6O1V 


1 + C, 


'^Pjn  _ 9Pjn  _ _ 3Pjn 
suns'  Silnil)  3ilnh 


= - (Po  - P, 


:r) 


1 + 0.601 


1 + 2Cf 


‘(f 


1 + C. 


'^erV'i'er^- 


(5.16) 


The  graphical  relation  between  pj^j^  and  S'  is  shown  in  Figure  5.4.  Eqs.(5.l6) 

would  be  used  in  the  analysis  of  externally  pressurized  gas  bearings  with 
inherent  restrictors. 


5.2  Centrally-Fed  Strip  (Figure  5.5) 

With  impermeable  surfaces  and  a uniform  gap.  Equation  (5.1)  reduces  to 
Integrating  repeatedly,  one  finds 


r- 1 -*d  ^ ' 
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3 ^2  1 


- liyijj 


(5.17) 


i(p2-p  2)  = 


12yi|;  ^ - 

3 (a-x) 


for  0 < a-x  < i 


ijj  is  the  pressure  flux  and  is  a constant  of  integration.  The  boundary 
condition  p(x=a)  = p^  has  been  imposed  in  obtaining  Equation  (5.17).  At 


the  edge  of  the  inherent  restrictor 


1 (p.  2 _ p 2)  = 

Eliminating  ip  between  Eqs.C5.l6)  and  (5.18),  one  obtains 


Vi2  ,,  24ufc^Ts 


■\/k(1  - •5^)1’  2 - K(1  - ^)1 
cry'  Ps  L Ps  J 


where  K 


= (1  _ Ecrfi  + 0 
i Ps  •- 


r(fijl)2  _ (^)2 

I^^Ps  ^ ^Ps^ 


.601V  1 + C, 


(5.18) 


(5.19) 


Eq.(5.19)  shows  that  for  a given  (Pj/Pg)’  how  (Pi^/Ps^  varies  with  h as 
illustrated  in  Figure  5.6. 


Pin<  Per 


(5.20) 


Sonic  flow  at  the  restrictor  throat  occurs  at  p^j,  - Per-  P®*"  Pin  < Per*  flow 

is  determined  entirely  by  the  sonic  throat.  Then,  instead  of  Eq.(5.19),  gap 
is  related  to  p^j^  by  setting  ij;  to  in  Eq.(5.18);  that  is 

j.-V“T^p"‘'2T  “ Pin  < Per  <5-?0) 

Beyond  the  throat,  prior  to  the  full  development  of  the  viscous  velocity  profile, 
there  would  be  an  extended  region  of  supersonic  over  expansion,  then  followed 
by  a complicated  "high  Mach  number  viscous  flow".  A full  account  of  the  pressure 
field  precisely  by  analytical  means  is  not  yet  possible.  For  this  reason, 
it  is  common  engineering  prudence  to  avoid  restrictor  choking. 
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If  one  desires  to  compute  the  stiffness  of  the  bearing,  a variational  pertur- 
bation with  respect  to  h would  be  performed  on  Eq.(5.18): 


Pin  ^Pin=  - 36^nh] 


(5.21) 


With  the  aid  of  Eq.(5.16),  6<lnii)  would  be  eliminated,  yielding 


6p. 


in  _ 


6h 


(5.22) 


2 P 


in 


L(P,  - P,p[l  + 0.601(i=L^}  2 (p.^2_  p^2)  J 


The  numerical  value  of  5pj^jj/6h  is  approximately  proportional  to  the  bearing 
stiffness.  As  shown  in  Figure  5.6,  it  may  be  maximized  for  given  (Pg,  Pg) 
through  a choice  of  h.  This  is  generally  regarded  as  the  ideal  design. 


If  p /p  is  large,  then  an  "exit  choking"  condition  can  develop.  This  occurs  when 

S 3 

Pg 

■1.19(^) 

Ps  for  Y = 


f (Pg/P^)  - 1 


or 


4' 


cr 


.1-35(3 

ps 


and  are  marked  in  Fig.  5.6  with  the  symbol  It  is  seen  that  the  optimum 

design  in  the  case  p,/p  = 5 is  not  valid  because  "exit  choking"  would  have 

S cl 

taken  place.  "Exit  choking"  can  be  avoided  by  discrete  feeding,  this  will  be 
discussed  in  Section  5.4. 


5.3  Various  Axially-Symmetric  Bearings 

The  above  development  can  ho  readily  generalized  for  various  other  bearing 
geometries  so  long  as  rotational  symmetry  exists  with  the  centerline  of 
the  supply  hole.  With  impermeable  bearing  walls,  referring  to  Chapter  3, 
the  governing  equations  for  common  bearing  configurations  are: 

Journal  Bearing 


fc  (h'p  = 0 


(5.23a) 
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Thrust  Bearing 

" fr  ^ 

Conical  Bearing 

sin>r  — (sin  rh^  ~ 

T r . I- 

Spheru-al  Bearing 

sini  (sint>iiV  -I-’)  " 

. !■ 


(5.23b) 


(5.23ci 


( 5 . 2 id ) 


If  the  film  thickness  is  unilnnn,  then  with  a siiitahlo  coordinate  transfor- 
mation each  of  tliese  four  equations  can  he  reduced  into  Equation  (5.17)  .so 
that  the  cor  resi^ond  ing  results — e.g.,  Eigure  5.6 — would  be  d i rec  t 1 c 
applicable,  provided  the  term  in  the  definititm  of  tlie  restriition 

parameters  be  accordingly  identified.  In  the  case  of  a journal  bearing, 
there  is  no  need  to  perform  coordinate  transformation;  one  mav  immediatelv 
recognize  the  axial  coordinate  "z"  to  be  equivalent  to  ”x"  in  Equation 
(5.15)  and  "f,"  would  simply  be  the  axial  distance  from  the  feeding  plane 
to  the  ambient  end.  For  thrus?  and  conical  bearings  Equation  (1.27)  provide-i 
tile  needed  transformat  ion,  while  F.(|uation  (3.28)  furnishes  the  I ransforn at  ion 
for  a spherical  bearing.  Thus,  referring  to  Figure  5.6,  the  equivalent 
for  these  bearings  arc  as  follows: 


Equivalent 


Journal 

E/2 

Thrust 

r . 1 itn(r  /r  . ) 

in'  a in 

Conical 

esc’’  r . 1 5n  (r  , 
in  a 

Spher  ical 

tanij)  72 
in 
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5.4  Feedins  Throu^;li  Discrfte  I lol  os 

In  the  preceding  sections,  flow  supply  is  assumed  to  enter  tiu  hearing  gap 
througti  a line-feed.  In  order  to  achieve  the  desired  degree  of  restriction 
without  "exit  choking"  (see  section  5.6),  the  "line-feed"  is  often  replaced 
hy  a series  of  small  feed-holes.  Although  qualitative  aspects  ot  the  flow 
process  through  the  discrc'te  feed-holes  are  essentially  unchangeil,  the  followin’ 
phenomena  may  sub.sl  an  t i.i  1 i y alter  the  restriction  behavior  quantitatively: 

• The  eftective  film  pressure,  p.  , at  each  feed-hole  may  he  peakid 

1211 

substantially  above  that  given  by  liquation  (5.18). 

• Upon  entering  the  bearing  I ilm,  flow  periphery  may  dilate  s ign i I ican ( I v 

I 22  i I 

prior  to  lull  di'Ve I opment  ol  the  viscous  velocity  profile. 

• Sonic  condition  is  more  like-ly  to  occur  at  the  throat  of  each  discrete 

restrictor.  In  the  event  ot  a sonic  inlet  ..  ondilion,  a supersonic 

over expans  ion  region  mav  cover  a large  enough  area  to  require  corri'ction 

1 2a] 

to  the  pressure  profile. 

If  tliere  are  n feed-holes  per  unit  Imigth  <'l  1 c’ed  line,  then  the  restrictor 
flow  area  is  '’tna  times  that  t>f  the  Mini  inlet  area,  "a"  being,  tf.c  tliroat 
radius  of  each  feed-hole.  If  n is  < boson  to  be  sufficiently  large,  pressure 
peaking  at  each  feed-hole  m.iv  be  mglected.  I'hen  the  analysis  ol  Seitions  5.2 
and  5.3  can  be  directly  used  for  the  discretely  fed  bearing  provided  V in 
eqs.  (5.11),  (5.13),  (5.14),  (5.15t  ind  (5.1b)  an.  rigarded  as  the  flux  through 
the  discrete  restr  ie  tt'rs . Thus,  in  eqs.  (5.17)  and  (5.18),  C shinild  hi’ 
multiplied  by  the  factor  (2nna).  ( oitsequetif  1 \ , eq.  (5.19)  i\tn  he  ad.jpled  by 

adding  a factor  (dtina)  to  its  right  hand  side.  The  giaplis  of  Mg.  5.b  would 
then  be  applicabli’  if  (Ztuia  ) ln'  divided  into  the  abscis.si  anil  mult  iplied  to 
(f/'Cer)’  symbol  on  each  curve  gives  the  upper  bound  of  (dnna)  to  avoid 

"exit  choking".  For  e.<ami>le,  for  a di.itomic  gas  with  rj.'Pj”5,  "exit  choking" 
is  avoided  by  requiring 

2Tina  < 0.2921 

When  this  is  done,  however,  the  restrictor  would  eiiienntor  sonic  ci’uiditioi' 
at  the  corresponding  point  so  that  eq.  (5.19)  .should  be  replaced  bv  eq.  (5..’0). 
Irrespective  of  the  question  of  su)iersonic  overexpans  i I'n  downstreams  ol  the 
restrictor  throat,  the  portions  for  H’  larger  than  that  marked  by  " + 
should  be  revised. 
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5.5  Non-uniform  Bearing  Gap 

In  the  event  of  a non-uniform  gap,  the  pressure  profile  would  be  al tired  in 
the  film,  even  if  the  boundary  pressures  (Pip,  p^)  'ire  unchanged.  This  lie- 
pendence  of  the  pressure  profile  on  the  bearing  gap  distribution  may  iiave 
a significant  influence  on  the  bearing  stiffness.  An  illustration  of  this 
issue  will  be  given  by  a detailed  analysis  of  a spherical  bearing,  whicii 
is  fed  through  a port  at  the  pole  witli  an  arbitrary  axial  eccentricity  as 

shown  in  Figure  5.7.  To  focus  our  attention  on  the  filir.  pressure  profile, 
restrictor  behavior  will  be  put  aside;  i.e.  p^^  (instead  of  Pg)  is  considered 
to  be  given. 


The  problem  at  hand  is  specified  by  the  following  set  of  equations: 


sin'<'h^p  ^ = -12uRt(; 
ci(|> 


h = C (1  - Ecosii) 


(5.24) 


'’('^in)  = Pin''  P(^>a>  = Pa 


This  problem  can  be  simplified  by  the  transformation 


n = 


/Cs3  d(t> 
sin* 


(5.25) 


y ielding 


dp  _ _ 12pRil; 


dn 


C3 


i(p2  . p^2,  . . (p  . 


hf.) 


(5.26) 


p2  _ p 2 P ~ ho 

p 2-p2  p _ri 

' in  a in  a 


where 


(i>)  = (i  -e^)  ^ 3f.^cosx  + cos^X  + e(3  + r/")  8n  sinx  + (1  + 3r^)  fn(tan^) 


x(if)  = 2 tan“  ( 


-1  / 1 + t 


1 - e 


tan 


(5.27) 


^ * * »■  ' . 
'mm  “ ■ 


^ d ^ 
- 


A 


1 
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Graphical  illustration  of  the  above  result  is  given  in  Figure  5.8.  It  is 
seen  that,  as  e increases  (gap  decreases),  the  area  under  the  pressure  profile 
decreases.  This  is  the  indication  of  a negative  stiffness!  Thus,  unless 
there  is  a sufficient  restrictor  compensation  of  this  effect,  this  type  of 
bearing  tends  to  "lock  up". 

The  tendency  for  "lock  up"  to  take  place  is  aggravated  by  large  eccentricity 
which  causes  significant  gap  non-uniformity.  It  is  more  serious  in  the  case 
of  a discretely-fed  bearing  than  that  of  a line-fed  bearing. 
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5.6  Exit  Choking 

For  a line-fed  strip  bearing,  sonic  flow  is  most  likely  to  orcnr  at  the  exit 
because  the  film  pressure  is  lowest  so  that  flow  speed  is  highest  tiiere.  In 
a rigorous  treatment  of  a liigh  "Mach  number"  gas  film,  one  should  make  allowance 
for  inertia  forces  in  the  momeuLum  balance  and  also  include  transverse  viscosity 
variation  due  to  temperature  variation.  In  order  to  gain  a first  insight  in 
the  overall  physical  situation,  the  problem  will  lirst  be  exaiuined  with  very 
rough  approximations  to  minimize  mathematic.il  drudgery.  Subsequently,  the 
same  problem  will  be  reviewed  with  .1  more  rigorous  treatment  and  t lie  basis 
for  a more  accurate  tlieoretical  treatment  will  he  accordingly  indii.iied. 

5.6.1  An  Approximate  Treatment 

For  the  present  purpose,  tiie  isoviscous  a pproxim.at  ion  will  be  ado[>ied,  ,ind 
the  convective  inertia  will  be  included  with  a presumed  paraliolic  velocity 
profile.  Such  a treatment  should  be  regarded  as  a first  appox  im.i  t ion  (see 
Sections  2.2  and  2.4). 


The  assumed  parabolic  velocity  profile  is 


Furthermore,  the  continuity  condition  is  approximated  as 


( 5.28) 


(5.20) 


Thus,  using  above  approximations  in  estimating  convective  inert  i.i,  then,  together 
with  the  perfect  gas  law  p = p/(RT),  one  finds  tiie  nu'mentum  b.il.'ince  to  yield 

^ 1 ri£  _ .2^; 

2 p L 3x  3x 

3y 

= iiE  [1  - lli  (5.  iO) 

P 3x  j^4 

The  second  term  on  the  right  hand  side  is  convective  inerti.i  effect  peculiar 
to  a compressible  fluid  since  it  originated  from  u ' . 
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Integrating  twice  with  the  boundary  conditions  u(0)  = u(li)  = 0,  one  obtains 

^ iLlZdll  AE  1 ^ (2y‘^-  4hy^+  h^y^+  h\-  + h"^) 

2w  3x  ^ 5AT  ,4  ' (5.31) 

h 


Returning  to  the  continuity  condition,  one  can  write 
h 


pudy  = const  = 


0 


aiT 


Or,  t'rom  eq.  (3.31 ) , 
3 


\2u 


ill  _ i iE  = 

3^  “ 55  p dx='  ^3 


(5.32) 


This  can  bt-  integrated  with  respect  to  x to  yield 


1-  2 *2^  *2,  ,p  , 

yip  -p  ) - P ^nr‘^) 


' 2jiL/  * ^ 

'—Y^(x  - X ) 

h 


where  the  superscript  "*"  refers  t('  the  choking  condition 
P '35  h.RT  a*h  a*-^  ' 


(5.33) 


(j.34) 


and 


“*  - f « >' 


(5.35) 


may  be  regarded  as  the  effective  sonic  speed.  Comparing  with  eq.(5.18),  one 
reckons  that  the  term  £n(p/p*)  in  eq.  (5.33)  represents  the  approximate 
contribution  of  convective  inertia. 

Using  symbols  defined  by  eqs.(5.34  and  (5.35),  eq.(5.32)  can  be  rewritten  as 


/£  _ J.  (£_)  = _ llHiL* 

p ^ dx  'p*-^  2 


(5.  ib) 


P*h 


For  any  p(x=0)  = eq.(5.36)  indicates  that  a moiU)ton ica  1 1 y d»-cri'asing 

pressure  distribution  cannot  be  continued  below  p* . Thus,  a correct  solutimi 
of  the  film  pressure  must  satisfy  the  Inequality 

p > p*  (5.  37) 
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Glven  the  half-wldtli  of  tlic  bearing  strip  i,  along  with  (u,  a*,  p.  , p ), 

X n ci 

eq.(5.36)  can  be  integrated  to  yield 


1 ,,^in^^  , ' ,p*,  „ ,^in^  12  a*l 

T ^ - 1 - (-^)  in  (- — ) = ^ 

2 P*  [Pa  ..  Pa  Pa  Pa’^ 


(5.38) 


which  may  be  solved  to  find  the  integration  constant  (p*/p  ),  which  must  not 

cl 

be  larger  than  unity  if  the  inequality  (5.37)  is  to  be  satisfied.  Violation 
of  the  latter  condition  can  occur  if  I is  too  small.  The  lower  bound  of  C 


is  found  by  setting  p*  = p^  in  eq.(5.38),  i.e. 


r h , p . T p . 

= '|-TT*  2 

I2pa*  z p p 

a a 


This  relationship  is  graphically  shown  in  Fig.  5.9.  If  i.  ' «,* , the  exit  film 

pressure  p*  satisfies  the  euijation 


1 2pa*  _ _1  ,P in 

p.  h2  ■ 2 
in 


P’'  P . 

( ) ?n  (^i) 

Pin  P" 


(5.40) 


The  apparent  discontinuity  between  p*  and  p^  is  a necessary  anomaly  of  tlie 
present  analysis. 


The  exit  clicking  phenomenon  is  seen  to  depend  on  two  parameters  which  may  be 
recognized  as 


The  lower  bound  of  the  former  parameter  is  fixed  by  eq.(5.39).  Kxlt  choking 
occurs  if 

• p /p.  is  too  small, 

*^a  '^in 

• I is  too  narrow,  or 

• h is  too  large. 

Representative  pressure  profiles,  with  and  without  exit  choking,  are  illustrated 

in  Fig.  5.10.  In  this  example,  exit  choking  takes  place  for  p.  /p  > 10. 

2 2 2 X n 3 

(P  “P  )/(p  "p  ) is  plotted  against  x/J..  In  the  absence  of  convective  inertia, 

3 S (X 

such  a plot  would  be  a straight  line.  With  p /p  = 6,  there  is  very  little 
evidence  of  any  significant  convective  Inertia  effect.  As  the  condition  of 
exit  choking  is  approached,  the  pressure  gradient  steepens.  Accordingly,  the 


-H  j- 


pressure  near  the  exit  would  increase  with  flow  increase  which  would  ct)rrcspuiul 
to  a larger  gap.  Thus,  it  is  clear  that  "exit  choking"  causes  a "negative 
stiffness"  tendancy. 


5.6.2.  A Rigorous  Re-examination 

Credibility  of  the  above  analysis  is  a bit  uncertain  since  the  exit  choking 

condition  in  fact  contradicts  the  apriori  assumption  of  an  isothermal  gas. 

[25] 

Fortunately,  a more  satisfactory  treatment  is  available  , which  is 

consistent  with  the  compressible  boundary  layer  theory  and  satisfies  the 
energy  equation  by  dint  of  the  relatively  inocuous  assumption  of  a unity 
Prandtl  number.  This  approach  will  be  briefly  outlined.  It  will  be  shown 
that  the  qualitative  aspects  of  the  approximate  analysis  presented  above 
are  not  contradicted,  albeit  the  accuracy  of  the  quantitative  results  may 
be  in  doubt. 


Adopting  the  boundary  layer  approximation  and  assuming  a unity  I’randtl 
number,  the  governing  equations  for  the  steady  state  flow  of  a gas  between 
parallel  plates  are 


lx  3y 


, 3 u , .1  u . 

(ur-  + V—)  = 
9x  3 V 


3 . 3u 

dt  ^ ^y) 


15.41) 


p = p<Rt 

xM  + li! 

Y-l  ^ 2 


= const. 


Y-l 


where  T^  is  the  wall  temperature  which  is  assumed  to  be  a constant. 


Now,  making  the  substitutions 
u = pu,  V = pv 


(5.42) 


one  f inds 


A A 

dU  3v 
3x  3v  " 


0 


(5.4  )> 
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''  9u  •'  3u 


- <1  - S ^ 4;  < Iv' " “ 


where 


^ [3  + 2(y  - 1)M^]  (y  - l)Ma 
[l  + (y  - DM^]^  YP 


(5.44) 


.2  = _iL 


. a2  = YiiT 


(5.45) 


In  this  form,  compressibility  effects  show  up  in  the  coefficient  of  dp/dx 
and  the  term  9.  Otherwise,  the  system  of  equations  is  exactly  same  as  that 
of  the  incompressible  entrance  flow.  Numerical  solution  reveals  that,  except 
in  the  entrance  region,  the  velocity  profile  is  approximately  a parabola 
even  when  the  Mach  number  is  quite  large. 


Eqs.  (5.43)  can  be  integrated  across  the  channel  to  yield 


ilE  / /-I  - ^ 

P dx  = 


v=h  -I 

[■  € - 

y=0  ( 


edy  - f- 
dx 


(5.46) 


It  may  be  noted  that  the  first  two  terms  on  the  right  hand  side  of  eq.(5.46) 
are  of  the  same  sign.  The  second  term  represents  the  influence  of  transverse 
density  variation  on  the  wall  shear  stresses  and  is  inversely  proportional 
to  the  film  pressure.  The  third  term  mainly  accounts  for  change  in  tlie  veloc- 
ity profile.  A detailed  numerical  study  of  eq.(5.43)  revealed  that,  beyond 
the  entrance  region,  the  velocity  profile  is  not  too  different  from  a parabola. 
One  can  thus  devise  a suitable  approximate  calculation  method  for  eq.(5.46) 

A 

by  the  assumption  of  a suitable  distribution  of  u.  Write 


u = u f (y/h) 
o 


(5.47) 


f is  the  distribution  function  such  that 


f d (f)  = 1 


(5.48) 


/ ' ' 

> «<■ 
r 


“ ^ ^ W ' ^ \ 
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u is  thus  the  mean  value  of  u.  One  can  then  define  a "Mean  Mach  Number"  to 
o 


O V Y p 


(5.49) 


Accordingly,  eq.  (5.46)  can  be  rewritten  as 


)h^  , 1 I 1 

o dp  1 df  I , 

nu~H'F  Pd^  = ~ 

w o**'  w I 

y=0  A 

^ h ,, 

+ 

U ax 
w 


where. 


T 

^o>  = f 

Jo 


^ (1  - M2)d(y/h) 


rOf  1 ° [3  + 2(y-1)m2]  [ df  1 y- 


(5.50) 


(5.51) 


(5.52) 


H = ^ f2d(y/h) 


(5.53) 


Upon  given  the  distribution  function  f and  the  viscosity-temperature  relation- 
ship u(T),  f’(Mj^)  and  G(M^)  can  be  calculated  as  universal  functions  of  M^. 

This  is  done  for 


(5.54) 


u = const. 

H is  found  to  be  a constant  and  is  precisely  1/10.  F and  G are  dependent  on 
both  Y and  as  shown  in  Figure  5.11.  F is  seen  to  vanish  at  M^  ' 1 for 
either  value  of  y.  In  fact,  F - 1 - M^^  would  be  a very  adequate  approximation 
for  the  diatomic  gas.  In  any  case,  the  vanishing  condition  of  F would  establish 
exit  choking.  From  this  point  of  view,  one  might  observe  that,  the  "isothermal 
approximation"  presented  previously  would  somewhat  prematurely  predict  choking  at 


/ ^ "TW*'  ' 
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which  is  0.68  for  y = 1.4.  On  the  otlier  hand,  C presents  an  additional 
effect  here  which  tends  to  mask  this  discrepancy  somewhat. 


A reasonable  approximation  for  F and  f!  Co  permit  closed  form  integration  of 
eq.  (5.50)  is 

* 2 


F = 1 - (--^) 
P 

* 

P-> 


(5. 55) 


* * 

where,  p ^ and  p ^ "’“’’y  take  on  the  following  numerical  values,  corresponding 
to  the  dashed  lines  shown  in  Figure  5.11: 


1.4 

1.6667 


/-X-  ' 

w o 
1.0 
0.95 


* 

I At  it 

W o 
0.8 
1.0 


where  upon,  eq.(5.50),  lor  a constant  h,  can  be  integrated  to  become 

(p,  - p)  [f  (p„  - P)  - p; ! - (p^ ' - p;  ^ in  ('^-i ‘’f 

P + P2 

12p  u P T X 
wo  w 


-) 


(5.56) 


where  p^  is  the  entrance  fluid  film  pressure  (neglecting  inherent  restriction) 


r\  * * 

at  X = 0.  Choking  takes  place  at  (p  = p^  , ,<  = ? );  i.  e., 
1 


^Po  ■ PP  (2  ^Po  ^ Pi  - P2  j 

. * 

12p  u T e 
wo  w 


A 

* 2 * 2 P + P, 

+ (?2  - Pj  ) fat 

Pi  + I’2 


h-' 


(5.57) 
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(p  = p , X = X.)  would  be  an  "unchoked"  solution  if  the  correspond  ins  u^, 

which  satisfies  eq.(5.56),  would  yield  i<  I from  eq.(5.57).  For  a 

* 

choked  exit  condition,  u is  found  from  eq.(5.57)  with  i = and  then 
p(o«:  X $ 1)  > p^.  These  features  are  qualitatively  unaltered  from  the 

previous  "isothermal  approximation". 


5 .7  Nomenclature 


(.'hapter  / 


Symbol 

Del  ini t ion 

Section 

a 

Half-width  of  inlet  slot 

5.2 

a 

Feed-hole  radius 

5.4 

a 

Acoustic  speed 

5.6.2 

* 

a 

Effective  sonic  speed 

5.6.1 

C 

o 

EinplricaJ  coefficient  in  eq.  (5.13) 

5.1 

F 

Function  ol  M , eq.  (5.51) 

5.6.2 

G 

Function  of  eq.  (5.52) 

5.6.2 

h 

Bearing  gap,  h^  - h^ 

K 

A constant , eq . (5.19) 

5.  1 

Width  of  strip  bearing  gap 

* 

i 

Lower  bound  of  i without  exit  choking 

M 

Mach  number 

M 

o 

Mean  Mach  number,  eq.  (5.49) 

5.6.2 

P 

Fluid  pressurt 

Pa 

Ambient  pressure 

Per 

Critical  pressure  for  isentropic  expansion 

5.1 

Peff 

Effective  film  inlet  pressure  of  incompressible  fluid 

5.1 

Pin 

Compressible  film  inlet  pressure 

Po 

Static  pressure  at  film  inlet 

5.1 

Ps 

Stagnation  pressure  of  supply 

Pth 

"Throat"  pressure  at  film  inlet 

5.1 

P* 

Pressure  at  choking 

5.6.1 

SI 

Gas  constant 

T 

Absolute  temperature 

c > 


Symbol 

T 

s 


T 

w 

u 

u* 

u 

u 


o 


U 

tn 


V 

V 

6W 

X 

X 

X 

y 

V 

5^ 


Def Inltlon 

Stagnation  temperature  of  supply 
"Throat"  temperature  at  inlet 
Temperature  of  bearing  wall 
x-component  of  velocity 

Velocity  profile  for  calculation  of  convective  inerti.i 
Mean  velocity 
pu 

Mean  value  of  u 

Mean  velocity 

Throat  velocity  at  inlet 

y-component  of  velocity 
pv 

Load  deficiency  due  to  entrance  pressure  drop 

Coordinate 

Distance  from  entry 

Pressure  drop  ratio,  eq.  (5.1A) 

Normal  coordinate  measured  from  reference  surface 
Ratio  of  specific  heats 

"Displacement  thickness"  of  entrance  pressure  profile 


2 Two-dimensional  surface  spatial  gradient  operator 

G Axial  eccentricity  ratio  of  spherical  bearing 

0 Compressible  dissipation  function,  eq.  (5.44) 

M Viscosity  coefficient 

u Value  of  u at  wall 

w 


Chapter/ 

Section 

5.  I 

5.1 

5.6.1 

5.6.1 

5.6.2 
5.6.2 


5.1 
5.  1 

5.  1 


5.6.2 

5.6.2 


V 


Kinematic  viscosity 
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Chapter/ 


Symbol 

Definition 

Section 

P 

Density  of  gas 

"th 

"Throat"  density  at  Inlet 

5.1 

4) 

Colatt itudinal  angle  of  spherical  bearing 

5.3,5.: 

4 

Inlet  mass  flux 

X 

"Sommerfeld  angle" 

5.5 

4' 

Pressure  flux 

T 

Dimensionless  pressure  flux,  eq.  (5.11) 

5. 1 

rr 

Value  of  \p  for  choking 

5.1 

6.0  TIME  DEPENDENT  EFFECTS 


Transient  variation  In  the  gas  bearing  film  pressure  is  a mass  storage 
effect  and  may  be  considered  as  compressibility  transient.  It  is  exhibited 
both  in  terms  of  an  initial  transient  and  also  in  terms  of  a "steady- 
state"  time-delay.  Mathematically,  as  shown  in  Eqiiati<'n  (2.24),  a first 
order  time-derivative  appears  simultaneously  with  second  order  space 
derivatives.  The  structure  of  this  equation  is  similar  to  the  diffusion 
equation  which  is  commonly  encountered  in  transient  heat  condition.  In 
principle,  the  problem  can  always  be  solved  by  a "marching"  process  numeri- 
cally. However,  this  typo  of  calculating  procedure  is  known  to  be  vulnera- 
ble to  numerical  instability;  therefore,  special  attention  must  be  given  to 
the  algorithm  in  order  to  realize  both  spatial  accuracy  and  temporal 
stability  [26]. 

Step-by-step  numerical  integration,  however,  is  not  always  a practically 
attractive  approach.  If  one  is  content  with  perturbation  solution  and  is 
not  concerned  with  the  initial  transient,  then  the  linearized  time  dependent 
problem  can  be  treated  in  the  frequency  domain.  In  Section  6.1  the  dynami- 
cally perturbed  gas  bearing  problem  will  he  rigorously  examined  in  the 
context  of  operational  calculus.  In  Section  6.2,  the  frequency  domain 
solution  will  be  examined  in  terms  of  its  global  characteristics. 

Because  compressibility  transient  is  basically  nonlinear,  there  is  a non- 
linear rectification  effect  when  large  amplitude  time  dependence  is  involved. 
Thus,  a gas  film  with  sustained  finite  amplitude  oscillation  is  capable  of 
supporting  steady-state  load.  The  theoretical  description  of  such  a 
phenomenon  is  treated  in  Section  6.3. 

6.1  Time-Dependent  Perturbation  Problems 

As  a model  problem,  the  time-dependent  perturbation  problem  will  be  examined 
in  some  detail  for  a self-acting  bearing  with  impermeable  surfaces.  For 
porous  bearings,  mass  storage  in  the  porous  material  must  be  also  accounted 
for.  Similarly,  for  externally-pressurized  bearings,  pressure  transient  in 
the  supply  line  must  be  included  if  upstream  flow  restrictions  are  present. 


r 
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The  governing  equations  were  previously  given  as  Equations  (2.27)  through 
(2.30).  To  facilitate  subsequent  discussion,  the  following  simplified 
notations  are  introduced  for  the  perturbation  problem: 

h ^ ^ 

} 6p  ^ {-  (U^  + U2)  '2/  (Pq«P)> 


.h  . Po-^Po  ^ 


p }Ah 
’ o 


(6.1) 


Then  Equation  (2.30)  can  be  rewritten  as 


-).5h 


(6.2) 


Let  the  boundary  conditions  be 


ip(r^,  t)  = 0 


(6.3) 


The  subscript  "i"  indicates  the  possibility  of  mure  tlin  one  bearing 
boundary.  I,et  the  initial  conditions  be 


5p(r,  t=0)  = n(r):  i^h(r,  t-0)  = H(r) 


(6.4) 


1271 

One  can  operate  on  Equation  (6.2)  with  the  I.aplaco  transformation^  so  th.it 


4p 


— S t. 

e qp  dt;  4li  = 


e ^ 6 h d t 


(6.5) 


Then  one  obtains 


(*+sh)4p  = hll+pH-(  +sp)6h 

o 0 0 o 


(6.6) 


and  the  boundary  conditions  become 


6p(r^)  = 0 


(6.7) 


In  principle,  one  may  regard  the  system  of  equations  governing  ftp  as 


L 


r ✓ * 


• 


r 


ar-  ^ ■-*  j ..  t 
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solvable  by  the  Green's  function  method.  Let  G(r;  u)  be  the  Green's 
function;  then  by  Its  definition. 


(>■  + h^S)G  = '(r,  li) 


(6.8) 


/(r,  u) 


Is  the  two-dimensional  Dirac  delta  function  defined  by 


0 if  r is  excluded  from  A 

u 

1 if  r is  included  In  A 

u 


(6.9) 


A^  Is  any  subdomain  of  the  bearing  surface.  u is  the  dummy  position  vector 
of  the  area  element  dA.  And 


G(r^,  u)  = 0 


(6.10) 


Accordingly, 


= G{h^n  + p^H  - ( + sp^)fhfdA  (h.Il) 

The  above  integral  covers  the  entire  bearing  surface.  All  t(-rms  enclosed 
in  brackets  are  regarded  as  tunctioiis  of  tlie  dur.jny  position  vi-ctor  u. 

Finally,  the  inverse  L.ip'aie  transfi'rm  c.in  be  "erformed  on  Enuation  (6.11) 
so  that 

6p(r,  t)  ’ '{4p(r,  s)} 


1 

2’'.i 


g-jon 


St 

e 


op  ds 


(6. 12) 


is  a suitable  constant  which  is  larger  than  the  real  part  of  any  polo  of 

<5p.  In  general.  Equation  (6.12)  indicates  the  sum  of  a number  of  exponential! 

tima  dependent  term.  The  exponential  coefficients  are  the  poles  .'f  6n.  Hut 

according  to  Equation  (6.11),  the  poles  of  6p  must  be  also  the  poles  of  Gfih 

since  the  spatial  integration  process  weighted  with  either  h n or  p H would 

o o 
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not  cause  any  singular  behavior.  If  the  system  is  stable,  then  the  initial 
conditions  associated  with  II  and  H would  decay  exponentially.  At  the  same 
time,  since  the  initial  conditions  do  not  affect  the  exponential  coeffi- 
cients, the  question  of  perturbed  dynamic  stability  is  Independent  of  the 
initial  conditions. 

For  a stable  system.  Equation  (6.12)  can  be  evaluated  with  6 ^ 0.  However, 
if  6h  contains  a purely  harmonic  time  dependence;  i.e.,  6h  - then  6h 

would  have  poles  at  s = + joj.  Tt  would  then  be  necessary  to  regard  the 
limiting  process  6 0 to  take  place  subsequent  to  the  performance  of  the 

Integra]  operation.  Equation  (6.12)  can  be  rewritten  as 

oo 

• ► 1 V > 

('■p(r,  t)  = — {f(r,  u;)cosu)t  - g(r,  uj)  s incut  ;dw  (6.13) 

^ o 

(f,  g)  are  respectively  the  real  and  imaginary  parts  of  6p(r,  s’-j'o).  If  p 
is  solved  with  both  (II,  H)  set  to  zero,  then  the  above  procedure  would  lead 
to  a special  solution  with  null  initial  conditions.  The  late  time  behavior 
of  any  solution  with  a different  set  of  initial  conditions  is  not  distinguish- 
able from  this  special  solution. 

6.2  Global  Characteristics  of  Harmc^nically-Perturhed  Gas  bearings 

Again,  tlie  self-acting  gas  bearing  with  Impermealile  walls  will  be  used  as  a 
model  problem.  Equations  (6.1),  (6.2),  and  (6.3)  would  be  the  system  of 
governing  equations.  in  the  use  of  operational  calculus  to  solve  a linear 
time-dependent  problem,  the  response  to  an  impulse  lunction  together  witii 
the  convolution  theorem  allows  one  to  compute  the  response  to  .jny  form  of  time 
dependent  driving  function  [27].  The  response  to  an  impulse  function  i an  he 
expressed  as  a Fourier  integral  of  harmonic  response  in  the  Inversion  formul.i. 
For  this  reason,  the  dynamic  environment  is  often  specified  in  the  frequency 
domain.  The  harmonic  perturbation  solution  is  also  a convenient  tool  in  the 
evaluation  of  dynamic  system  stability  in  accordance  with  an  energy  flciw 
criteria  [28,29].  Margin  of  dynamic  stability  is  one  of  the  most  important 
issues  in  the  design  of  high-speed  rotating  equipment.  Thus,  the  harmonic 
perturbation  analysis  of  a gas  bearing  is  a subject  of  considerable  impor t inco . 
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The  liarinonic  perturbation  solution  can  in  principle  be  express<>d  in  terms 
of  a Green's  function.  For  computation  purposes,  the  Green's  function 
formulation  for  gas  bearings,  however,  cannot  be  conveniently  implemented 
due  to  the  complexity  in  the  two-dimensional  operator  V'  as  given  by 
Equation  (6.1).  Therefore,  the  hannonic  perturbation  problem  will  be  recast 
in  a form  which  can  be  ri'adily  solved  by  relatively  conventional  numerical 
procedures . 

To  specify  a harmonic  perturbation  problem,  the  perturbed  be.iring  gap  is 
specified  by 

6h(r,  t)  = 6g(r)  cosvt 


6g(r)  describes  the  spatial  distribution  ol  the  pi'r turbat  ion , and  v is  the 
particular  frequency  under  consideration.  The  cosine  I unction  is  chosen 
to  indicate  that  time  referenci*  is  fixed  by  the  gap  motion.  With  the 
omission  of  initial  conditions,  the  dynamic  perturbation  (irohleni  is  homo- 
geneous in  (6p,  5h) . Consequently,  6p  must  also  have  a h.irmonic  time- 
dependence  at  the  same  frequency  v.  Furthermore,  one  can  write 

ilh  = 6g(r)  Re{e^'’''^l  (6.11) 

and  the  unknown  perturbed  pressure  may  be  expressed  as 

6p  = Ref  oq(r)e^'^*' } (6.1(0 

i5q(r)  would  generally  be  complex,  indicating  the.  possible  presence  of  phast' 
shift  between  6p  and  6h. 

As  both  bearing  surfaces  are  assumed  to  he  fixed  to  rigid  bodies,  fig  may 
represent  the  effect  of  either  a translation  or  a rotation  of  one  surface 
relative  to  the  other.  in  the  case  of  a translational  perturbation,  the 
relative  motion  between  the  two  bodies  is  indicated  by  the  perturbation 
vector  fieRefe^^*'},  then  the  perturbation  of  tlie  bearing  gap  is 

= n • fie  Refe  } 


fib 


((1.  171 
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Thus,  by  comparison  with  Equation  (6.15), 

i5g  = n • '5c  (6.18) 

In  the  case  oC  un  angular  or  a rotational  perturbation  represented  by  the 
relative  angular  perturbation  vector  , bearing  gap  perturbation 

becomes 

6h  = - n • r X (5y  Re{e^^*'}  (6.19) 

Accord ingly , 

6g  = - n • r X oY  = r X n • 6Y  (6.20) 

With  6g  represented  by  either  Equation  (6.18)  or  (6.20),  substitution  ol 
Equation  (6.15)  and  subsequent  elimination  of  the  common  factors  e^^*^  and 
Ke{  } I rom  all  terms,  Eqiiat  ii'ns  (6 . I ) through  (6.3)  i>ecome 

- 1 - ► ^ 

1 > - 

= - .2  • ^2  0-^  + (-I-P^^ll^g)  - ivp^^dg  .n.21) 

<5q(r^)  = 0 ((1.22) 


Well-proven  numerical  methods  for  solving  tin'  above  system  of  cqui  t i on;-,  ir 
available.  Formally,  one  may  write 

6q  = ■'  '5g  (6.21) 

with  ••  indicating  a linear  operator.  Then  the  corresponding  per  t iirba  1 ion  of 
the  force  vector  is,  according  to  Equation  (2.33) 

= n 0>q  dA  = n(  6g)  dA  (b.?4) 
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For  the  translational  perturbation  problem,  making  use  of  Ffiuation 
one  can  write 


6F 


n (:'  n)  dA 


(Sc 


(6.25) 


Thus,  staying  with  the  vector  lUJtation,  one  can  define  the  I ransl at ion.i  1 
stiffness  matrix  to  be 


force 


n ( n)  dA 


(6.2b) 


For  the  most  general  case,  a perturbation  force  may  also  be  induced  bv  an 
angular  perturbation.  However,  for  most  commonly  used  he.irings,  there  is 
a sufficient  degree  of  symmetry  so  that  it  is  usually  possib|<.  to  locate 
the  origin  of  the  coordinate  system  to  eliminate  indu(.  tion  of  fori  e by 

angular  perturbation  as  well  as  induction  of  moment  bv  translational  j 

perturbation.  Therefore,  perturbation  of  the  moment  vector  will  next  be  j 

examined  with  Equation  (6.20),  yielding  i 


n (,. 


n)  dA 


6y 


(6.27) 


Again,  the  angular  stiffness  matrix  can  be  defined  as 


moment 


r X n (^r  x n)  dA 


(6.28) 


As  presented  in  Equations  (6.26)  and  (6.28),  there  are  no  apparent  formal 
distinction  between  a gas  bearing  and  an  incompressible  fluid-film  hearing. 
There  is  actually  a very  important  difference  which  is  implicit  in  the 
structure  of  the  operator  c,  which  reflects  the  lefthand  side  of  Equation 
(6.21).  For  an  incompressible  fluid-film  bearing,  they  are  all  related  to 
the  spatial  differential  operators  and  are  all  real.  The  righthand  side 
contains  a frequency  Independent  real  part  and  an  Imaginary  part  which  is 
directly  proportional  to  the  frequency.  For  this  reason,  the  "dynamic 
stiffness"  matrix  of  an  incompressible  fluid-film  bearing  is  usually 
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represented  by  a set  ol  "stiffness  iiuitrix  coefficients"  and  <i  set  of  "damping 
matrix  coefficients."  For  a gas  bearing,  because  of  the  presence  of  iv1i^6q 
on  the  lefthand  side  of  Equation  (6.2J),  both  real  and  imaginary  parts  of 
are  frequency  dependent.  Consequently,  one  can  regard  both  "stiffness  and  damping 
matrix  coefficients"  of  a gas  bearing  to  be  frequency  dependent. 

6.3  Rectification  Effects  of  Higli-Ft  equency  Oscillation.-. 

A gas  film  undergoing  sustained  oscillation  develops  a mean  pressure  above 
ambient.  Historically,  this  phenomenon  was  reported  as  evidence  of  an 
anomalous  rheological  property  of  air  [30].  Never tlieless , a correct  theore- 
tical analysis  as  well  as  the  recognition  of  the  true  meclianism  actually 
predated  the  intense  gas  lubrication  efforts  of  tlie  60's  [11]. 

By  sustaining  tlie  liigh-f requency  oscillation  witii  an  electrical  or  a magnetic 
means,  load  support  capability  can  be  developed  by  a tliin  gas  film.  A gas 
bearing  operating  on  this  principle  is  called  a squeeze-film  gas  bearing. 
Squeeze-film  gas  bearings  of  various  geometrical  shapes  are  illustrated  in 
Figure  6.1.  The  as>’mplf>tic  tiieory  for  squeeze-film  gas  bearings  is  briefly 
outlined  below. 

Consider  the  most  general  equation  of  the  isothermal  gas  bearing  equation 
for  impermeable  walls: 

3 ^ 

^ (ph)  f 2 (Wj  + P'2  P}  = 0 (6.29) 

If  a sustained  fluctuation  at  a sufficiently  high  frequencv  is  present,  the 
leading  term  of  Equation  (6.29)  dominates:  thus 

^ (ph)  ~ 0 (6.10) 

o C 

Time  independence  of  (ph),  however,  is  obviously  not  valid  at  bearing  edges; 
therefore.  Equation  (6.30)  is  a legitimate  approximation  only  in  the  internal 
region.  Suppose 


h 


0 


(6.31) 
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Various  Squeeze-Film  Gas  Bearings 
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Then  Equation  (6.29)  is  reduced  to 


24y 


^ 2 = 0 


(6.32) 


If  one  performs  a "time-series  average"  on  Equation  (6.32),  one  obtains 


2 3 

■ 2 . < p h > = 0 


(6.33) 


The  symbol  < > designates  the  averaged  result  of  the  enclosed  quantity. 

Since  p and  h are  both  uniform  along  the  bearing  edges. 


2 3 2 3 

< p h > = p <h  > 
^ a 


(6.34) 


is  a solution  of  Equation  (6.33).  This  result  may  be  regarded  is  an 
analogy  of  the  "mass  content  rule"  discussed  in  Section  4.4.2.  In  the 
internal  region,  then 

2 3 2 2 3 

■ p~h  > = (ph)“  < h>  = p^^  < li  > 


. < h^^ 

Ph  = Pa  -hT 


(6.33) 


If  h is  a simple  harmonic  function;  e.g.,  h = C(1  - ecosvt),  then 


'3  2 

ph  = 1 + 2 e P^P: 


(6.36) 


The  temporally  averaged  pressure  is  then 


3 2 

<p>  = 1+:te  p 


a l+Ecosvt" 


1 + I 


(6.37) 


If  E 1,  a significant  amount  of  pressure  elevation  above  p is  seen 

a 

possible.  The  question  remains,  how  high  should  the  frequency  be  to 
justify  Equation  (6.30)?  The  answer  lies  in  rewriting  Equation  (6.29) 
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In  a dimensionless  form,  scalingv^wlth  1/R,  h with  C,  and  p with  pg ; one 
would  then  establish 


12uvR 

a = — >>  1 


(6.38) 


PaC 


to  be  the  required  condition.  To  be  thoroughly  rigorous,  one  must  be 
cautious  that  Equation  (6.38)  does  not  violate  the  prerequisite  for  an 
isothermal  analysis;  namely. 


yC 


(-|^)(Y-1)M  « 1 


(6.39) 


Above  results  have  been  general ized  to  relax  the  assumptions  of  Equations 
(6.31)  and  (6.38)  [32,33],  and  to  consider  "long  time  scale"  dynamics  [34] 

as  well  as  "hybrid  squeeze-slide"  gas  films  [35). 
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6. A Nomenclature 


Symbol 

C 


->■ 

6e 


f 

6F 

S 

G 

li 

n 

5h 

6h 

H 


= force 


= monen  t 

M 

6M 

P 

-> 

n 


o 

6p 

6p 

? 

5q 

-1? 

-> 

r 

► 

r. 


Definition 


Ciiauter  / 
Sect  ion 


Specific  bent  at  constant  pressure 
Translaticuial  perturbation  vector 

A 

Real  part  of  6p  for  s-^ju) 

Perturbation  force  vector 
Ini,i4;inary  part  of  for  S'-jjj 

Spatial  clistr  ibut  ii'n  of  harmonic  gap  perturbation 

Green's  1 unction  of  i5p 

St e.idy-sta to  film  thickness 

Gap  piTiurl'.it  ion 

Laplace  transform  of  Ah 

Ini t la  1 c end i 1 1 on  of  i h 

Force-stiffness  matrix 

Moment-stiffness  matrix 

Linear  operator 

Mach  number 

Perturbation  moment  vector  due  to  6p 

Outward  unit  normal  vector  of  bearing  surface 

Ambient  pressure 

Steady-state  pressure 

Pressure  perturbation 

Laplace  transform  of  6p 

Linear  operator  for  fp,  eq.  (6.1) 

Spatial  distribution  of  harmonic  pressure 
Linear  operator  for  6p,  eq.  (6.1) 

Surface  position  vector 
Position  vector  of  ith  boundary 
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Chapter/ 

Def Ini 1 1 on  Sec  t ion 

Laplace  parameter 
t ime 

• > 

Dummy  variable  of  r 

Two-dimensional  Dirac  delta  function 
Angular  perturbation  vector 

Two-dimensional  surface  spatial  gradient  operator 
Excursion  ratio  of  harmonic  motion 
Thermal  conductivity 
Viscosity  coefficient 
Frequency  of  harmonic  motion 
Initial  condition  of  6p 
12pvR^ 

Squeeze  number,  

P 

a 


Subscripts 

1 referring  to  lower  bearing  surface 

2 referring  to  upper  bearing  surface 

<>  Time-series  average 
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